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ABSTRACT 

Experimental  data  for  the  charging  of  an  air  receiver  is 
presented  and  interpreted  in  detail.   The  data  indicates  a 
substantial  departure  from  the  adiabatic  behavior.   The  ex- 
perimental results  are  used  to  evaluate  existing  closed  form 
expressions  for  the  thermodynamic  state  of  a  gas  in  a  re- 
ceiver.  A  method  for  experimentally  determining  the  convec- 
tive  heat  transfer  coefficient  is  developed,  evaluated  and 
used  in  conjunction  with  these  expressions. 

The  experimental  work  was  performed  from  March  19  69 
through  May  19  69  at  the  Naval  Postgraduate  School,  Monterey, 
California. 
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I.   INTRODUCTION 

Systems  using  the  charging  of  a  receiver  are  often  em- 
ployed today.   A  common  practice  when  analyzing  a  receiver 
being  charged  is  to  assume  that  the  process  is  adiabatic. 
This  simplifies  the  governing  equations  for  the  system  and 
leads  to  a  fairly  straightforward  expression  for  the  state  of 
the  gas  in  the  vessel.   During  the  course  of  the  charging 
process  however,  a  substantial  temperature  difference  between 
the  walls  of  the  receiver  and  the  gas  in  the  receiver  may 
develop.   Often  the  tank's  thermal  capacitance  is  of  suffi- 
cient magnitude  to  permit  the  extraction  of  large  amounts  of 
energy  from  the  gas  with  only  a  small  change  in  the  tempera- 
ture at  the  outside  surface  of  the  tank  walls.   This  obscures 
the  apparent  effect  of  heat  transfer,  for  the  absence  of  a 
noticeable  temperature  change  at  the  outside  of  the  receiver 
walls  may  lend  credence  to  the  adiabatic  assumption  whereas 
actually  large  amounts  of  heat  are  being  transferred  from  the 
gas  to  the  receiver  walls.   Errors  introduced  by  such  an 
assumption  may  lead  to  serious  problems.   For  example,  if 
heat  trans fe]  Ls  neglected,  calculations  for  determining  the 
amount  of  gas  needed  to  charge  a  vessel  iav  lead  to  a  pre- 
dicted value  lower  than  that  actually  required. 

In  order  to  aid  in  the  analysis  of  this  problem  Reynolds 
[Ref.  1]  developed  a  theory  which  includes  the  effects  of 
heat  transfer  for  determining  the  thermodynamic  state  of  gas 
in  a  receiver  during  the  charging  process.   He  developed  four 
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closed  form  solutions  for  the  various  magnitudes  of  the  system 
parameters.  Using  his  criteria  to  determine  the  proper  closed 
form  solution,  many  charging  processes  may  be  analyzed  with- 
out excessive  difficulty.  The  accuracy  of  these  solutions  de- 
pends to  a  great  extent  on  the  accuracy  with  which  the  average 
convective  heat  transfer  coefficients  can  be  estimated  for  the 
charging  process. 

In  his  experimental  investigation  of  the  blowdown  process 
[Ref.  2]  Reynolds  tested  his  theory  using  an  h  (average  con- 
vective heat  transfer  coefficient  for  the  entire  blowdown 
process)  determined  by  taking  the  value  of  his  heat  transfer 
parameter  NTU  (number  of  thermal  units,  see  Appendix  A)  which 
when  used  in  his  theoretical  equation  of  state  gave  the  best 
fit  to  the  experimental  data.   These  h  values  were  on  the 
order  of  those  predicted  on  the  basis  of  an  assumption  of 
steady  state  turbulent  free  convection  inside  the  receiver, 
i.e., 

£|  =  .13(GRPR)1/3 

Therefore  Reynolds  suggested  the  use  of  an  h  based  on  this 
assumption  in  his  solution  to  the  charging  process. 

The  problem  still  remains,  however,  that  unless  h  is  de- 
termined by  an  independent  means  it  is  difficult  to  come  to 
any  quantitative  conclusions  as  to  the  accuracy  of  Reynolds1 
theory. 

The  purpose  of  this  study  has  been  to  derive  a  method  for 
experimentally  determining  the  average  convective  heat  trans- 
fer coefficient  between  the  gas  and  the  receiver,  and  to 
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apply  these  results  to  three  of  Reynolds '  closed  form  solu- 
tions in  order  to  evaluate  their  accuracy.   The  three  closed 
form  solutions  evaluated  are  for  conditions  where  the  heat 
transfer  from  the  receiver  walls  to  the  surrounding  medium 
can  be  neglected.   These  three  cases  are  the  ones  most  com- 
monly encountered.   The  fourth  of  these  closed  form  solutions 
deals  with  the  case  where  the  walls  of  the  receiver  are  very 
thin  and  the  heat  transfer  from  the  receiver  walls  to  the  am- 
bient medium  must  be  taken  into  account.   This  solution  was 
not  evaluated  due  to  the  difficulties  in  experimentally  deter- 
mining the  average  convective  heat  transfer  coefficient  be- 
tween the  walls  of  the  receiver  and  the  surrounding  medium. 
Also  presented  in  this  paper  are  representative  values  for 
the  various  heat  transfer  parameters  determined  in  this 
undertaking. 

In  the  following  sections  of  this  thesis  the  theoretical 
and  analytical  methods  and  results  are  discussed.   Following 
a  description  of  the  charging  system,  the  theoretical  model 
due  to  Reynolds  is  reviewed  and  the  methods  for  the  experi- 
mental evaluation  of  the  heat  transfer  coefficient  are  devel- 
oped.  The  experimental  results  are  then  discussed  in  the 
light  of  their  relevance  to  Reynolds'  model  and  their  useful- 
ness in  engineering  problems.   A  final  section  is  included  in 
order  to  summarize  this  work  and  indicate  some  areas  for 
future  investigation. 
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II.   OBJECTIVES 

Experiments  on  the  charging  of  a  gas  receiver  were  con- 
ducted at  varying  flow  rates  and  heat  transfer  environments, 
in  order  to  accomplish  the  following  objectives: 

(1)  Evaluate  the  closed  form  solutions  developed  by 
Reynolds  to  approximate  the  thermodynamic  state 
of  a  gas  during  charging. 

(2)  Develop,  use,  and  evaluate  an  expression  for  ex- 
perimentally determining  the  averaqe  instantaneous 
heat  transfer  coefficient  in  a  gas  receiver  being 
charged. 

III.   DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

A.   GENERAL  DESCRIPTION 

The  test  apparatus  (Figures  1  through  4)  consisted  of  an 
aluminum  cylindrical  gas  receiver  with  a  volume  of  1.03  cubic 
feet.   The  complete  physical  dimensions  of  the  test  apparatus 
are  listed  in  Table  1.   The  vessel  was  designed  according  to 
the  ASME  Boiler  and  Pressure  Vessel  Code  for  a  working  pres- 
sure of  200  psig.   Its  top  was  removable  to  allow  access  to 
the  interior  of  the  tank,  using  an  O-Ring  flange  arrangement 
to  insure  airtightness  when  the  top  was  in  place.   The  re- 
ceiver was  fitted  with  a  flow  metering  device  mounted  in  the 
center  of  the  tank  top.   The  other  fittings  consisted  of  a 
mounting  for  a  thermocouple  probe,  a  pressure  transducer 
mounting,  and  a  bleed  valve  arrangement. 
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This  test  receiver  was  mounted  so  as  to  be  immersed  in  a 
larger  tank  which  was  either  open  to  the  atmosphere  or  filled 
with  ice  water  depending  on  the  particular  run.   Air  for 
charging  the  test  tank  was  supplied  from  two  117  cubic  foot 
air  vessels  charged  by  an  air-cooled  compressor.   The  supply 
pressure  to  the  flow  metering  device  was  varied  from  180  psiq 
to  10  psig. 

B.   FLOW  METERING 

The  mass  flow  rate  for  each  run  was  determined  by  maintain- 
ing a  critical  flow  through  one  of  two  nozzles  or  throuqh  an 
orifice.   These  small  diameter  (D)  devices  were  machined  from 
stainless  steel  in  accordance  with  a  paper  by  Grace  and  Lapple 
[Ref.  3]  and  the  discharge  coefficients  (Cn)  given  in  their 
paper  were  used  in  these  calculations.   The  diameter  and  dis- 
charge coefficients  are  listed  in  Table  1.   The  orifice  or 
nozzle  was  mounted  on  the  tank  top  in  a  flange  type  arrange- 
ment that  was  sealed  by  O-Rings  and  connected  to  a  1  inch  line 
that  lead  to  the  supply  tanks.   The  stagnation  pressure  (P~) 
was  read  by  a  local  pressure  tap  and  gage  arrangement  mounted 
just  upstream  of  the  flow  metering  device.   The  stacrnation 
temperature  (TQ)  was  likewise  found  by  using  a  thermocouple 
mounted  in  the  flow  just  upstream  of  the  metering  device.   It 
should  be  noted  here  that  for  velocities  of  the  magnitude  oc- 
curring in  these  tests,  it  was  safe  to  assume  that  the  differ- 
ence between  stagnation  and  local  conditions  upstream  of  the 
orifice  was  negligible.   Thus,  assuming  adiabatic  flow  throucrh 
the  orifice  or  nozzle,  the  temperature  and  pressure  values  of 
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P_  and  TQ  indicated  by  the  apparatus  just  described  can  be 
used  in  the  critical  flow  equation, 

w  =  .532  C^  ~ 

The  air  to  the  flow  metering  device  was  controlled  by  a  cruick 
acting  Jamesbury  ball  valve. 

C.   TEMPERATURE  MEASUREMENT 

The  test  vessel  was  fitted  with  a  thermocouple  probe  con- 
sisting of  four  40-gage  copper  constantan  thermocouples.   The 
thermocouple  housings  were  designed  so  as  to  measure  as 
closely  as  possible  the  local  temperature  at  four  equally 
spaced  heights  in  the  receiver  (Fig.  1) .   The  thermocouples 
were  arranged  so  that  the  four  could  be  read  in  series,  and, 
at  the  same  time,  any  one  of  the  four  could  be  read  independ- 
ently.  The  series  signal  was  recorded  on  channel  one  of  a 
two  channel  Moseley  strip  chart  recorder,  while  the  signal  of 
the  thermocouple  being  read  singly  was  recorded  on  a  contin- 
uous Brown  recorder.   According  to  the  manufacturer's  speci- 
fications thermocouples  of  this  type  are  accurate  to  within 
±1.5°R  for  the  temperature  ranges  found  in  these  tests. 
Since  the  thermocouple  probe  readings  were  only  used  qualita- 
tively, no  further  calibration  was  made.   As  was  mentioned 
above,  the  stagnation  temperature  of  the  inlet  air  was  read 
from  a  single  thermocouple.   The  signal  of  this  40-gage  copper 
constantan  thermocouple  was  monitored  on  a  Leeds  and  Northrup 
recorder. 
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D.   PRESSURE  MEASUREMENTS 

The  pressure  in  the  test  receiver  was  obtained  by  use  of  a 
Daystrom  Wiancko  pressure  transducer  attached  to  the  receiver 
whose  signal  was  recorded  on  channel  two  of  the  Moseley  strip 
chart  recorder.   The  input  voltage  (21  volts)  to  the  trans- 
ducer was  supplied  by  a  Philco  power  supply  and  the  transducer 
output  signal  was  adjusted  so  as  to  read  10  psi/in  on  the 
Moseley  recorder.   The  pressure  transducer  was  calibrated 
originally  using  a  test  gage  and  later  using  a  dead-weight 
tester.   When  connected  to  the  recorder  the  pressure  trans- 
ducer registered  pressures  within  ±  1  psi  of  the  actual  pres- 
sure throughout  the  range  of  these  tests. 

IV.   EXPERIMENTAL  PROCEDURE 

The  experimental  program  was  divided  into  the  following 
three  cases: 

1.   In  the  first  group  of  runs  the  receiver  was  charged 
12  times  under  conditions  corresponding  to  the  ad- 
iabatic  case  in  Reynolds'  theory.   The  runs  con- 
sisted of  three  charges  at  each  of  four  possible 
supply  tank  pressures;  190,  180,  170  and  160  psig. 
The  3/8  inch  knife-edge  orifice  way  mounted  for 
these  tests  and  the  total  elapsed  time  for  a  run 
was  approximately  one  second.   The  large  tank  en- 
circling the  test  receiver  was  open  to  the  atmos- 
phere.  Thus  the  rapid  chargings  gave  NTU  values 
which  dictated  the  use  of  Reynolds '  adiabatic 
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model  for  theoretically  determining  the  thermodynamic 
state  of  the  gas  in  the  receiver  throughout  the  run 
(See  Section  V-A-3 ) . 

2.  The  second  group  of  runs  corresponded  to  the  isother- 
mal case  described  in  Reynolds'  theory.   The  36  runs 
consisted  of  3  runs  at  each  of  12  possible  supply 
tank  pressures  ranging  from  180  psig  to  10  psig,  in 
various  intervals.   The  1/32  inch  nozzle  was  mounted 
in  the  flow  metering  apparatus  and  the  run  time  var- 
ied from  150  seconds  to  75  seconds.   Once  again  the 
large  tank  surrounding  the  test  tank  was  open  to  the 
atmosphere.   The  resulting  very  small  mass  flow  rates 
were  introduced  to  correspond  to  NTU  values  in  the 
region  of  7,  a  lower  limit  for  Reynolds'  criteria. 

3.  The  last  set  of  runs  was  made  with  the  test  receiver 
completely  immersed  in  ice  water  contained  by  the 
large  tank.   The  first  9  runs  were  made  with  the 
1/32  inch  nozzle  in  place,  three  each  at  suoply  pres- 
sures of  180,  160  and  140  psig.   The  last  18  runs  of 
this  case  were  made  with  the  1/8  inch  nozzle  mounted 
in  the  tank  top,  three  each  at  supply  pressures  of 
180,  150,  120,  90,  60  and  30  psig.   The  ice  bath  was 
to  maintain  the  receiver  walls  at  a  constant  tempera- 
ture and  thus  present  the  gas  with  an  isothermal  sink 
to  correspond  to  the  one  described  in  the  analytical 
model. 
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The  values  of  the  supply  pressure  listed  above  reoresent 
the  nominal  charging  pressure.   The  actual  staanation  pressure 
just  upstream  of  the  flow  metering  device  was  recorded  for 
each  run. 

In  making  a  series  of  runs ,  the  supply  tanks  were  first 
charged  to  approximately  19  0  psig.   The  system  was  then  iso- 
lated from  the  compressor  and  the  test  receiver  was  allowed 
to  reach  equilibrium  with  its  surrounding  medium.   The  quick 
acting  valve  between  the  supply  tanks  and  the  flow  metering 
device  was  then  opened.   It  was  kept  open  until  the  pressure 
in  the  test  tank  reached  a  value  insufficient  to  maintain  a 
critical  pressure  ratio  across  the  flow  meter.   The  valve  was 
then  closed  and  a  run  terminated.   For  the  succeeding  run  the 
supply  tank  pressure  was  adjusted  to  the  desired  value  and 
the  test  receiver  was  vented  to  the  atmosphere.   The  next  run 
did  not  commence  until  the  temperature  of  the  gas  in  the  re- 
ceiver (as  indicated  by  the  thermocouple  probe)  reached  equi- 
librium with  that  of  the  tank  walls  and  thus  the  surrounding 
medium.   Having  made  one  pass  through  the  supply  pressures 
indicated  for  a  particular  case,  the  supply  tanks  were  re- 
charged and  the  procedure  repeated  for  a  second  and  third 
time. 

Throughout  the  charging  process  the  pressure  and  tempera- 
ture history  of  the  gas  in  the  tank  was  recorded  on  a  two 
channel  strip  chart  recorder. 
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V.   THEORETICAL  CONSIDERATIONS 

1 
A.   THEORY  DEVELOPED  BY  REYNOLDS 

1.   Introduction 

In  an  analysis  of  the  charging  process  one  object  is 
to  obtain  the  time  dependent  thermodynamic  state  of  the  eras 
in  the  receiver.   A  common  practice  is  to  express  the  tempera- 
ture of  the  gas  (the  dependent  variable)  as  a  function  of  the 
mass  (the  independent  variable) .   For  a  constant  volume  re- 
ceiver the  perfect  gas  law  may  then  be  used  to  express  the 
pressure  as  a  function  of  the  mass.   Knowing  the  mass  flow 
rate  the  thermodynamic  state  of  the  gas  as  a  function  of  time 
can  be  described. 

With  this  in  mind  Reynolds  develooed  the  following 
model  for  a  receiver  being  charged: 


///////// 


77777777/ 


The  thermal  capacitance  of  the  receiver  walls  and  anv  other 
internal  structure  is  lumped  into  a  single  caoacitance 


The  nomenclature  used  in  this  section  is  taken  from  Ref- 
erence 1,  see  Appendix  A  for  listing. 
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represented  by  C   and  equal  to  the  product  of  the  mass  of  the 

capacitance  and  its  specific  heat  (C   =  M  c  ) .   The  heat  trans- 
*  r  c     c  p 

fer  resistance  between  the  receiver  walls  and  the  surrounding 
environment  is  represented  by  R^  and  the  heat  transfer  resist- 
ance between  the  capacitant  material  and  the  qas  in  the  re- 
ceiver is  represented  by  R.. 

2.   General  Differential  Equation  for  Charging 

Reynolds  assumes  that  the  thermal  resistance  as  well 
as  the  thermal  capacitance  are  invariant  with  time.   He  also 
assumes  perfect  mixing  of  the  injected  gas  and  the  gas  in  the 
receiver,  and  that  the  walls  of  the  receiver  are  at  a  uniform 
temperature  throughout.   Using  these  assumptions  he  extends 
his  model  to  include  all  the  significant  enerqv  terms  as 
follows : 


////////A 


4     M,T 

/ 


I,  7        m 

^l^ide 


dU 

de 


///////// 

S777 


Combining  an  energy  balance  for  the  above  model  and  the  heat 
transfer  rate  equations, 

T-T 


^oo  = 


T  -T 
c   °° 

R 


<*i  = 


R. 


the  following  non-dimensional  general  differential  equation 
for  the  charging  is  obtained  (see  Appendix  B) : 
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d2T*  dw*     NTU+NTU^ 

w*M*  ^-4r  +  w*[2w*  +  NTU  +  M*  §Sr  +   — =-* M*]§sr 

dM*  dM*       CQ*        dM* 

,  *     ,  NTU  NTU 

+  [W*  d¥*"  +  c^(NTU+NTUoo)w*  +      t   °°]T* 

-  *    NTU+NTU  NTU  NTU 

-  kT*w*  ^St ?r-i — -  w*kT* =-, T*  =  0  . 

1    dM*      CQ*         1     C  *       °° 

For  the  case  of  constant  mass  flow  rate  (w*  =  1)  this  equation 
reduces  to: 

d2T*  dT* 

M*  2_J_  +  [C  +C0M*]5ut  +  C,T*  +  C,  =  0        (1) 
dM*       1   2    dM*     3       4 

where  C.  ,  C~  ,  C^,  and  C.  are  constants  comprised  of  the  sys- 
tem parameters  NTU,  NTU  and  C  *;  defined  as, 

x       (hA) 

NTU  =  =-= = , 

R.  c  wA     w  c 
l  v  0      o  v 

,       (hA) 

00 

NTU   = 


R  c  wA    wAc.T 
00  v  0     0  v 

and 


M_c.     C 
C 


*  -   c  P  ..  c 


0     Mnc     MAC 
0  v     0  v 

The  parameter  NTU  is  used  by  Reynolds  to  represent  the  conduct- 
ance (reciprocal  of  resistance)  of  the  receiver  as  well  as  to 
give  a  measure  of  the  rate  of  the  process.   The  CQ*  parameter 
is  the  ratio  of  the  capacitance  of  the  receiver  to  the  initial 
capacitance  of  the  gas. 

Reynolds  solves  equation  (1)  analytically,  however 
the  resulting  series  solution  is  seen  to  be  impractical  for 
many  engineering  applications.   Fortunately  in  manv  cases  the 
values  of  the  system  parameters  are  such  that  a  simpler  closed 
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form  solution  may  be  obtained.   These  solutions  can  be  obtain- 
ed by  simplification  of  the  series  solution  or  by  returning 
to  the  general  differential  equation  for  charging.   In  his  re- 
port Reynolds  nonetheless  chooses  to  develop  each  closed  form 
solution  from  a  simplified  form  of  his  original  model.   A 
brief  description  of  these  and  their  solutions  are  presented 
in  the  following  four  sections.   A  complete  derivation  of  the 
general  differential  equation  for  charging  as  well  as  deriva- 
tions of  the  closed  form  solutions  are  presented  in  Appendix 
B. 

3.   Adiabatic  Charging 

As  has  been  mentioned,  a  common  method  of  analysis  of 
a  system  being  charged  is  to  assume  the  process  is  adiabatic. 
This  assumption  of  no  heat  transfer  to  or  from  the  gas  may  be 
found  to  represent  the  case  of  a  high  mass  flow  rate  charging 
(wQ  large)  process  fairly  accurately.   This  may  be  thought  of 
as  due  to  the  fact  that  the  gas  has  not  had  time  to  transfer 
a  substantial  amount  of  heat  to  the  walls  before  the  charging 
process  is  terminated.   Even  if  the  high  mass  flow  rate  proc- 
ess is  fairly  lengthy,  the  effect  of  the  heat  transfer  is 
small  compared  to  the  effect  of  introducing  large  quantities 
of  mass  during  the  charge.   Such  an  assumption  may  also  be 
justified  for  the  case  where  the  gas  has  been  insulated  from 
the  thermal  capacitance  of  the  receiver  walls  (h  -  0) .   This 
case  corresponds  to  the  situation  where  NTU  is  verv  small, 
and  the  following  model  applies : 
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///A 


h    dM 


/////////A 

A 


dU  =  d(Mu) 


M,T 


The  solutions  for  the  dimensionless  temperature  and  pressure 

are , 

kT*-l 

T*  =  kTi  "  -w—  (2) 

and 

P*  =  kT*(M*-l)  +  1   .  (3) 

4 .   Isothermal  Charging 

A  receiver  with  a  large  thermal  capacitance  (C 
large)  and  either  charged  very  slowly  (w*  small)  or  having  a 
high  hA  value,  will  exhibit  an  isothermal  behaviour.  The 
normal  temperature  increase  in  the  gas  due  to  comoression  is 
not  observed  because  the  energy  is  extracted  from  the  gas  bv 
the  thermal  capacitance  of  the  receiver  before  it  can  become 
significant.  As  described  above,  a  system  of  this  type  will 
have  large  values  for  the  system  parameters  NTU  and  CQ*,  and 
may  be  represented  by  the  following  model: 


////////- 
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The  temperature  ratio  for  charging  is  by  definition, 

T*  =  1 
and  therefore, 

P*  =  M* 
5.   Charging  at  Constant  Mass  Flow  with  Heat  Transfer 
to  an  Isothermal  Sink 

In  some  systems  the  thermal  capacitance  of  the  wall 
is  much  larger  than  the  thermal  capacitance  of  the  eras  (C~* 
very  large).   For  this  case  it  can  reasonably  be  assumed  that 
the  energy  extracted  from  the  gas  is  not  great  enough  to  sig- 
nificantly affect  the  temperature  of  the  receiver  walls  (T 
is  a  constant) .   Thus  the  analytical  model  is  modified  to 
describe  the  heat  transfer  between  the  gas  and  an  isothermal 
sink  as  follows : 


///// 


/////////M 
777777777? 
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The  solutions  for  this  case  are, 


T*  = 


kT*+NTUT*  -  (kT*-l-NTU+NTUT*)M* 
1    *  c       1 c 

1+NTU 


-(1+NTU) 


(4) 


and 


P*  = 


(kT*+NTUT*)M*  -  (kT*-l-NTU+NTUT*)M* 
1 c_         1  c' 

1+NTU 


-NTU- 


(5) 


25 


6.   Charging  at  Constant  Mass  Flow  with  Inside  Resistance 

Negligible 

A  system  having  a  high  conductance  can  be  assumed  to 
have  negligible  inside  heat  transfer  resistance  (NTU  high) . 
This  implies  that  the  temperature  of  the  capacitance  is  the 
same  as  that  of  the  gas  and  the  heat  transfer  is  between  the 
capacitance  and  the  surrounding  environment.   That  is  distin- 
guished from  the  isothermal  case  in  that  here  the  thermal 
capacitance  is  finite  but  not  large  (very  thin  walled  cylinder 
for  example) .   Since  systems  of  this  nature  are  not  too  com- 
mon, Reynolds  feels  that  the  greatest  value  of  this  solution 
is  that  it  supplies  information  for  determining  the  effect  of 
capacitance  on  the  heat  transfer  in  a  charging  process.   The 
model  for  this  system  is: 


/     dU         d(Mu) 

„„„„A  de  =  de 

7  7        3M 

/'/  ///////  A 

/ 
/ 
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The  solutions  for  this  case  are 


1+NTU 


T*  = 


C*+l 
(1+NTU  -kT*-NTU  T*)\C*+M* 

OO  00   00   \   0 

1+NTU 


kT*  +  NTU  T* 

1  00   00 


(6) 
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and 


P*  =  i" 


1+NTU  -kT*-NTU  T* 

OO  00   oo 

(C0+1)  Wl+NTU  )T*-kT*-NTU  T* 

oo'  oo  ■ 


_1 3-1 

1+NTU 


-C* 


(7) 


B.   CRITERIA  FOR  APPLICATION  OF  CLOSED  FORM  SOLUTIONS 

In  the  previous  four  sections,  closed  form  solutions  were 
developed  for  the  state  of  the  gas  in  a  receiver  being  charged, 
These  solutions  were  developed  by  simolifving  the  original 
model  with  assumptions  as  to  the  magnitude  and  imoortance  of 
various  terms  in  the  general  solution.   For  example,  in  the 
adiabatic  charging  case  the  parameter  NTU  is  assumed  to  be 
very  small,  thus  the  mass  flow  rate  is  large  comoared  to  the 
heat  transfer  coefficient  and  a  solution  based  on  no  heat 
transfer  is  formulated.   Similar  tvpes  of  assumptions  are 
made  for  the  other  three  cases.   Reynolds  then  produces  quan- 
titative criteria  for  the  use  of  such  assumptions.   These 
criteria  are  based  on  a  maximum  deviation  of  5%  in  T*  at  the 
value  of  M*  =  9.   They  are  determined  by  comparing  the  case 
in  question  with  the  case  which  would  give  the  maximum  depar- 
ture from  this  behavior.   Thus  the  adiabatic  case,  (no  heat 
transfer  to  the  walls)  is  compared  to  the  case  where  the  ca- 
pacitance of  the  receiver  walls  is  infinite,  the  isothermal 
sink  case  (large  amounts  of  heat  transfer  to  the  walls). 
When  the  percentage  T*  difference  between  the  two  solutions 
is  plotted  against  NTU,  the  value  of  the  two  solutions  agree 
within  5%  for  values  of  NTU  less  than  .25,   In  a  similar 


z  / 


fashion  criteria  are  developed  for  all  the  closed  form  solu- 
tions.  Some  solutions  have  more  than  one  possible  criteria 
for  their  use,  however,  only  the  criteria  for  the  three  cases 
tested  here  which  are  met  by  our  exoerimental  system  are 
listed  below.   As  mentioned  before,  the  three  closed  form 
solutions  evaluated  here  are  the  most  useful  and  the  criteria 
for  their  application  are  those  most  commonlv  found. 
Adiabatic:   0  <  NTU  <  0.25,  all  C* 
Isothermal  Sink:   0.25  <  NTU  <  7    C*  >  40 
Isothermal:   NTU  >  7    C*  >  40 

C.   DETERMINATION  OF  CONVECTIVE  HEAT  TRANSFER  COEFFICIENT 

In  order  to  test  Reynolds'  theory  on  the  state  of  a  gas 
in  a  receiver  being  charged,  a  technique  was  devised  to  obtain 
a  value  for  the  convective  heat  transfer  coefficient  in  the 
receiver  in  which  the  tests  were  run.   To  this  end  the  follow- 
ing model  of  the  experimental  system  was  constructed: 


H.m 

l 
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where , 

m  =  the  mass  flow  rate  into  the  receiver 
q  =  the  heat  transfer  to  the  receiver  walls 
d(Mu)/dt  =  the  rate  of  change  of  the  internal  energy 
of  the  gas  in  the  receiver  with  respect 
to  time 
H.  =  the  stagnation  enthalpy  of  the  entering  gas 
T  =  the  temperature  of  the  gas  in  the  receiver 
P  =  the  pressure  of  the  gas  in  the  receiver 
M  =  the  mass  of  the  gas  in  the  receiver 
Several  points  should  be  noted  here.   First  it  was  assumed 
that  the  temperature  throughout  the  receiver  could  be  repre- 
sented by  a  single  value  T  (See  discussion  of  imperfect  mixing 
section  V-E) .   The  same  was  true  for  the  pressure  term  P.   It 
should  also  be  noted  that  this  development  does  not  include 
the  case  where  there  is  heat  transfer  between  the  surrounding 
medium  and  the  walls  of  the  receiver. 

An  energy  balance  on  the  receiver  yields, 

d(Mu) 
mH .  -  q  =    ... 
1    ^      dt 

The  kinetic  energy  of  the  gas  in  the  receiver  was  neglected 
for  it  can  be  shown  that  for  the  flow  rates  in  this  experi- 
ment, the  velocity  of  the  gas  in  the  receiver  was  verv  small 
when  compared  to  the  internal  energy.  Using  the  notation  f 
to  represent  the  derivative  of  a  function  f  with  respect  to 
time,  the  above  equation  can  be  written  as, 


ihH .  -  q  =  if 


itiu  +  Mu 
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Assuming  the  gas  in  the  receiver  is  thermallv  and  calorif- 

ically  perfect, 

q  =  -m(c  T-c  T.)  -Mc  T  (8) 

^        v    p  1       v 

where  T.  is  the  inlet  stagnation  temperature  of  the  gas.   By 

definition  of  the  convective  heat  transfer  coefficient, 

hA(T-T  )  =  q  =  -MC  T  -  ih(c  T-c  T.) 
w     M       v        v    D  1 

In  the  same  manner  as  T  represents  the  average  gas  temperature, 

T   is  used  to  represent  the  average  temperature  of  the  re- 
w 

ceiver  walls.   Thus  the  h  defined  here  renresents  an  average 
instantaneous  heat  transfer  coefficient  for  all  points  on  the 
inside  of  the  receiver  walls.   The  symbol  A  represents  the 
total  inside  area  of  the  receiver  walls.   By  use  of  the  per- 
fect gas  law  it  can  be  seen  that, 

p  .  m   T 

P    M    T 

for  a  constant  volume  system.   Therefore  the  above  equation 
can  be  written  as , 

-Mc  T(^-)+mc  T. 

h  =  v;  P   ■  P  x 

A  T-T 
w 

or,    using   the    perfect   gas    law   again, 

-c   P  (S+mc   T. 

The  next  step  was  to  determine  an  expression  for  T   if  it 

L  w 

was  not  maintained  a  constant  as  in  the  isothermal  sink  case. 
In  doing  so  the  difficult  task  of  determining  exDerimentally 
an  accurate  average  wall  temperature  was  avoided.   Since,  for 
simplicity,  no  heat  transfer  from  the  outside  of  the  receiver 
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walls  to  the  surrounding  medium  was  assumed,  an  enerqy  balance 
on  the  walls  yields, 

Q  =  4r(M  u  )  =  M  CT 
^    dt   w  w     ww 

where  C  is  the  specific  heat  of  the  walls  and  is  assumed  con- 
stant.  Combining  this  with  equation  (8) , 

q  =  -m(c  T-c  T.)  -  Mc  T  =  M  CT 
^        v    p  1       v      ww 

is  obtained.   Inteqrating  this  equation  over  time  from  t  =  0 

to  t  =  t  gives , 

/*;=A  M  CT  dt  =  |^=^[-m(c  T-c  T . ) -Mc  T]dt 
J t=0   w   w      Jt=0      v    D  1     v 

or 

M  C(T  -T   )  =  /*;=A  -3r(Mc  T)dt  +  fjrj:  mc  T.dt   . 
ww  w'    J  t=0   dt   v       j  t=0   d  i 

In  these  tests  T.  was  a  constant,  therefore, 

i 

M  C(T  -T   )  =  -Mc  T  +  MAc  TA  +  (M-MA)c  T. 
w    w   wA        v      0  v  0        0   D  1 

where  subscript  0  denotes  the  condition  at  time  t  =  0. 

Solving  for  T   gives 

3       w  * 

-Mc  T+MAc  TA+Mc  T,-Mnc  T. 

T   =  2 °  V  Q    P  X       S  P  I    +  T 

W  MwC  w0 

Again  using  the  perfect  gas  law  to  eliminate  the  term  T,  the 

final  form  of  the  expression  for  T   is  obtained, 

w 

T»  =  T„0  +  fc)  [|:VT1  -  fire)  [p-V    •    (10> 

0    \  w  /   0         \  w  ' 
Equations  (9)  and  (10)  furnish  a  means  for  estimating  the 
film  coefficient  h  in  terms  of  the  tank  pressure  and  gas  mass 
and  their  rates.   The  massive  experimental  sirmolicity  thus 
introduced  is  the  major  justification  for  the  acceptance  of 
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inaccuracies  stemming  from  the  definition  of  the  heat  transfer 
in  terms  of  spacial  averages  of  gas  and  wall  proDerties. 

D.   DISCUSSION  OF  ASSUMPTION  OF  NO  EXTERNAL  HEAT  TRANSFER 

In  the  derivation  of  the  expression  for  h  it  was  assumed 
that  there  was  no  heat  transfer  from  the  outside  of  the  re- 
ceiver walls  to  the  surrounding  medium.   This  simplifies  the 
energy  balance,  for  if  there  is  no  heat  transfer  from  the  walls 
to  the  outside  medium,  all  the  energy  leaving  the  gas  must  be 
stored  in  the  receiver  walls.   This  assumption  is  not  onlv 
valid  for  this  experimental  system,  but  for  many  pressure  ves- 
sels used  in  charging  processes.   A  vessel  capable  of  with- 
standing a  substantial  pressure  is  likely  to  be  constructed 
with  materials  and  dimensions  that  give  it  a  fairly  large 
thermal  capacitance  when  compared  to  that  of  the  entering 
gas.   In  the  test  receiver,  for  example,  the  thermal  capac- 
itance of  the  receiver  was  approximatelv  8.7  Btu/°R  while 
even  at  its  highest  mass  the  gas  capacitance  was  only  of  the 
order  of  .08  Btu/°R.   Thus  the  gas  temperature  drop  due  to 
heat  transfer  to  the  receiver  walls  would  have  to  be  extremely 
large  before  any  significant  change  in  the  wall  temperature 
was  observed. 

This  assumption  of  no  external  heat  transfer  may  also  be 
justified  by  considering  the  length  of  time  it  would  take  the 
temperature  increase  due  to  heat  transfer  at  the  interior  sur- 
face to  reach  the  outside  walls.   One  method  of  analyzing  this 
transient  conduction  heat  transfer  problem  would  be  to  consider 
the  receiver  walls  a  semi-infinite  solid  with  heat  transfer 
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due  to  convection  at  the  surface.   Thus  by  determining  the 
time  (t)  required  for  an  observable  temperature  change  to 
occur  in  the  semi-infinite  solid  at  a  distance  (x)  from  the 
surface  equal  to  the  receiver  wall  thickness  of  3/8  inch,  some 
idea  of  the  temperature  change  at  the  outside  surface  of  the 
test  tank  will  be  obtained.   According  to  Car slaw  and  Jaeger 
[Ref.  4],  the  solution  for  this  type  of  problem  is  given  by 
the  equation, 


v  _hx+h2 at    _„.«,_/  X         .    ,   it- H 

V 


—  =   erfc -e  erfc +   hvat\ 

\2/at/  12/at  / 


where 

a  =   the  thermal  diffusivity  of  the  material 

h   =  h/k  =  the  convective  heat  transfer  coefficient/ 

the  thermal  conductivity 
v  =  the  difference  between  the  wall  temperature  at 
time  t  and  a  distance  x  from  the  surface  and  the 
initial  wall  temperature 
V  =  the  difference  between  the  gas  temperature  and  the 
initial  wall  temperature. 
Using  the  highest  value  of  V  recorded  throughout  our  runs 
(100°R)  and  a  corresponding  value  for  h  we  find  that  v  is 
still  only  approximately  1.5°R  even  after  2  minutes  have 
elapsed.   Thus  for  our  runs  the  assumption  of  no  external  heat 
transfer  during  a  charge  is  justified.   This  will  be  the  case 
for  all  receivers  of  thermal  capacitance  sufficient  to  store 
the  thermal  energy,  transferred  from  the  gas,  with  negligible 
increase  in  temperature. 
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E.   DISCUSSION  OF  THE  EFFECTS  OF  IMPERFECT  MIXING 

In  the  derivation  for  determining  h,  and  in  the  closed 
form  solutions  derived  by  Reynolds,  it  was  assumed  that  the 
effects  of  imperfect  mixing  could  be  neglected.   In  other 
words,  it  was  assumed  that  for  analytical  purnoses  the  gas 
can  be  considered  perfectly  mixed  so  that  no  temoerature  or 
pressure  gradients  exist.   As  was  seen  in  an  unpublished  ex- 
periment on  the  heat  transfer  in  a  closed  container  (with 
similar  dimensions  to  the  test  receiver)  after  gas  injection 
[Ref.  5] ,  this  is  not  actually  the  case.   Temperature  gra- 
dients do  exist  and  this  leads  to  uncertainty  in  evaluating 
the  heat  transfer  data  correctly.   In  order  to  determine  the 
effective  temperature  potential  for  heat  transfer,  use  must 
be  made  of  some  sort  of  average  gas  temperature  so  as  to  keep 
the  complexity  of  the  experimental  equipment  and  the  data  re- 
duction process  reasonable.   If  this  average  gas  temperature 
is  dependent  upon  the  degree  of  mixing,  then  the  assumption 
that  the  gas  is  perfectly  mixed  may  lead  to  erroneous 
conclusions . 

In  his  paper  [Ref.  1] ,  Reynolds  showed  through  the  use  of 
a  simple  comparison  between  mixed  and  ummixed  gas  systems  that 
the  pressure  is  independent  of  the  degree  of  mixing  in  an  ad- 
iabatic  receiver  regardless  of  its  magnitude.   The  magnitude 
of  heat  transfer  in  these  tests  was  never  very  large  and, 
since  the  average  gas  temperature  used  in  the  derivation  for 
determining  h  is  based  directly  on  the  perfect  gas  law,  the 
temperature  thus  calculated  is  fairly  independent  of  the  de- 
gree of  mixing. 
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Reynolds  also  showed  that  the  temoerature  averaqed  with 
respect  to  mass  in  such  an  adiabatic  receiver  was  independent 
of  the  degree  of  mixing.   Realizing  that  this  average  temper- 
ature would  be  very  difficult  to  find  experimentally,  he  went 
on  to  compare  the  temperature  averaged  with  respect  to  mass 
and  that  averaged  with  respect  to  volume  and  found  that  the 
discrepancy  was  so  small  "that  use  of  a  volume  average  tem- 
perature in  experimental  investigation  is  entirely  satis fac- 

2 
tory. "    Therefore,  the  average  temperature  measured  by  the 

four  thermocouples  connected  in  series  at  equal  volumes 
vertically  in  the  tank  could  be  used  as  the  temperature  refer- 
red to  as  the  gas  temperature  in  the  closed  form  solutions. 

VI.   SUMMARY  OF  EXPERIMENTAL  RESULTS 

A.   FORM  OF  RESULTS 

The  data  from  a  particular  run  consisted  of  a  pressure  and 
temperature  history  for  the  test  receiver  during  charging,  and 
values  for  the  inlet  stagnation  pressure  and  temperature. 
From  this  data,  using  the  critical  flow  equation,  the  mass 
flow  rate  for  each  run  was  calculated  and  thus  the  mass  of  the 
gas  could  be  determined  for  any  instant  of  time.   Using  this 
mass  and  the  pressure  history  of  the  run,  the  instantaneous 
h  values  were  determined  by  equations  (9)  and  (10) .   The  val- 
ues of  NTU  were  then  calculated  and  used  to  determine  which 
class  of  closed  form  solution  the  criteria  would  designate 


2 

Reference  1,  p.  85 
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for  this  run.   The  values  of  P*  =  P/Pn  as  determined  bv  the 
appropriate  closed  form  solution  were  then  compared  with  the 
experimental  values.   The  average  percentage  deviation  be- 
tween these  two  values  and  the  average  value  for  NTU  through- 
out a  run  were  then  tabulated.   Plots  were  made  comoarinq  the 
P*  values  as  predicted  by  Reynolds  and  those  observed  experi- 
mentally for  representative  runs.   These  results  are  in  the 
form  of  P*  versus  M*  =  M/M-.  because  the  pressure  readings 
were  slightly  more  accurate  due  to  the  uncertainties  involved 
in  determining  the  temperature  of  a  gas  in  motion. 

The  values  of  h  used  xn  determining  NTU  have  been  repre- 
sented by  plots  of  the  Nusselt  number  (NU  =  hL/kf)  versus 
the  product  of  the  Grashof  and  Prandtl  (GRPR  =  L3p^.g3fAT/ 
yi(c  y/k)  _)  numbers  for  various  runs.   The  subscript  f  here 
indicates  that  the  term  was  evaluated  at  the  film  temoerature, 
which  for  this  case  was  the  mean  temperature  between  that  of 
the  gas  and  the  wall.   The  characteristic  length  L  of  the 
system  was  taken  to  be  the  receiver  height  for  these 
calculations . 

B.   HIGH  MASS  FLOW  RATE  RUNS 

A  series  of  runs  was  conducted  with  the  test  receiver  ex- 
posed to  the  ambient  atmosphere  at  high  mass  flow  rates  so  as 
to  correspond  to  the  adiabatic  closed  form  solution  develooed 
by  Reynolds.   The  dimensionless  capacitance  term 

C*  =  C  /M-.C  ,  used  in  the  analytical  criteria  was  approx- 
0    c'  0  v  J 

imately  657  for  all  runs.   Based  on  the  criteria  MTU  <  .25 

for  any  value  of  C*,  all  of  these  runs  were  found  to  correspond 
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to  the  so  called  "adiabatic"  case.   Values  of  1.4  for  the 
ratio  of  specific  heats  of  the  gas  (air  in  this  case)  and  1.00 
for  the  non-dimensional  inlet  stagnation  term  T* ,  (the  inlet 
stagnation  temperature  of  the  gas ,  T,  /  the  original  tempera- 
ture of  the  gas,  Tn)  were  used  in  calculating  the  theoretical 
P*  values  from  equation  (3).   Figure  5  is  a  graphical  compar- 
ison of  the  actual  and  theoretical  P*  values  for  representa- 
tive runs  of  this  type.   NU  versus  GRPR  plots  for  a  few 
typical  runs  of  this  nature  are  seen  in  figures  11,  12,  13 
and  14 ,  the  points  having  been  plotted  at  1/10  second  inter- 
vals throughout  the  runs.   Table  2  summarizes  the  results  of 
this  series  of  tests. 

C.   LOW  MASS  FLOW  RATE  RUNS 

A  second  group  of  runs  was  conducted  at  low  mass  flow 
rates  so  as  to  correspond  to  the  isothermal  closed  form  solu- 
tion proposed  by  Reynolds.   Once  again,  the  test  receiver  was 
exposed  to  the  surrounding  atmospheric  conditions.   Reynolds1 
criteria  for  an  isothermal  charge  solution  to  hold  is  NTU  >  7 
and  C*  >  40.   The  value  of  C*  for  this  group  of  runs  was 
again  approximately  657  and  the  NTU  values  calculated  varied 
from  1.8  to  14.6.   Thus  these  runs  straddle  the  cutoff  value 
of  NTU  =  7  and  can  be  used  to  determine  the  validity  of  this 
figure.   The  solution  for  an  isothermal  charge  is  simply 
p*  =  m*,  for  T  is  a  constant  and  therefore  T*  =  1.   The  val- 
ues of  P*  predicted  by  Reynolds  and  those  observed  experi- 
mentally are  compared  graphically  in  figure  6  for  represent- 
ative runs  of  this  type.   The  h  values  found  in  a  few 
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typical  runs  are  presented  by  the  NU  versus  GRPR  graphs  of 
figure  15,  16,  17  and  18.   The  points  in  these  figures  repre- 
sent 15  second  intervals  throughout  the  runs.   Table  3  sum- 
marizes the  results  of  this  series  of  tests. 

D.   ISOTHERMAL  SINK  RUNS 

The  final  group  of  runs  was  made,  with  the  test  receiver 
immersed  in  an  ice  bath,  at  intermediate  mass  flow  rates  so 
as  to  correspond  to  Reynolds'  isothermal  sink  model.   The 
criteria  for  this  case  is  simply  C*  >  40.   However,  in  order 
to  avoid  crossing  into  the  adiabatic  or  isothermal  solution 
regions,  the  value  of  NTU  must  be  greater  than  .25  and  less 
than  7  respectively.   It  should  be  understood  that  the  ice 
bath  was  not  necessary  in  order  to  meet  Reynolds'  criteria: 
with  or  without  it  the  value  of  C*  (approximately  629  for  this 
case)  was  large  enough  to  dictate  the  isothermal  sink  solution, 
The  ice  bath  was  employed  to  ensure  that  the  temperature  of 
the  receiver  walls  was  a  constant  and  that  a  significant  tem- 
perature potential  existed  between  the  receiver  walls  and  the 
gas.   Values  of  1.4  for  the  ratio  of  the  specific  heats  of  the 
gas  and  1.00  for  the  dimensionless  inlet  temperature,  T* ,  as 
well  as  for  the  dimensionless  wall  temperature,  T*,  (the  re- 
ceiver  wall  temperature  T  /  the  original  temperature  of  the 
gas  T„)  were  used  in  calculating  the  theoretical  P*  values 
from  equation  (5) .   Figure  7  is  a  graphical  comparison  of  the 
actual  and  theoretical  P*  values  for  representative  runs  of 
this  type.   NU  versus  GRPR  plots  for  a  few  typical  runs  of 
this  nature  are  seen  in  figure  19,  20  and  21,   The  points 
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were  determined  at  15  second  intervals  for  figure  19,  and  1 
second  intervals  for  the  other  two  figures.   Table  4  sum- 
marizes the  results  of  this  series  of  tests. 

VII.   DISCUSSION  OF  RESULTS 

In  order  to  systematically  analyze  the  test  results,  the 
following  section  is  divided  into  three  separate  sections. 
The  first  section  is  devoted  to  an  analysis  of  the  equation 
derived  for  the  convective  heat  transfer  coefficient.   The 
results  obtained  when  comparing  the  experimentally  observed 
values  of  P*  to  those  predicted  by  Reynolds  are  then  discus- 
sed.  The  last  section  deals  with  an  analysis  of  the  convec- 
tive heat  transfer  coefficients  obtained  with  emphasis  on  the 
identification  of  general  trends  that  might  be  investigated 
in  a  more  comprehensive  study  on  this  topic  alone. 
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A.   ANALYSIS  OF  THE  EXPRESSION  FOR  DETERMINING  THE  CONVECTIVE 
HEAT  TRANSFER  COEFFICIENT  EXPERIMENTALLY 
In  an  earlier  section  the  exoression 

-c  P  (=-)+mc  T. 

h  =  pT7 (9> 

a  r  (~)  -t  l 

was  derived  with  the  value  for  T  ,  when  it  was  not  considered 

w 

a  constant,  given  by  equation  (10). 

Tw  -  T„n  +  te)  'snTi-Tl  -  £rc)  tp-po]        (10) 

0    \   w  /   0         \   w  ' 
As  was  mentioned  before,  these  expressions  were  purooselv 
developed  to  obtain  a  value  for  h  from  experimentally  deter- 
mined data.   Due  to  the  difficulty  of  measuring  the  effective 
average  temperature  of  the  gas ,  this  term  was  eliminated  in 
the  above  expressions.   The  equation  for  h  is  thus  a  function 
of  P  and  M,  the  other  terms  being  system  parameters  which  were 
generally  held  constant  and  known  to  a  high  degree  of  accuracv 
throughout  a  run. 

First,  considering  the  expression  for  the  wall  temperature, 
equation  (10) ,  it  is  seen  that  even  when  the  temperature  of 
the  receiver  was  not  held  constant  the  temperature  of  the  walls 
changed  very  little  during  a  run.   This  can  be  seen  by  noting 
the  magnitude  of  the  terms  in  parentheses  as  well  as  the  fact 
that  the  pressure  and  mass  of  the  gas  are  never  more  than  one 
order  of  magnitude  larger  than  their  original  conditions  (MQ , 

P0) .   An  investigation  of  the  values  calculated  for  T   shows 
u  w 

that  this  was  indeed  the  case  for  these  tests. 
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Next  the  expression  for  the  heat  transfer  coefficient  it- 
self will  be  considered.   Since  the  derivation  of  this  exores- 
sion  was  centered  about  the  fact  that, 

h  =  A(T-T  J 
w 

it  can  be  seen  that  the  numerator  of  expression  (9)  is  an  ex- 
pression for  the  heat  transfer  q  and  the  denominator  is  the 
product  of  the  inside  area  of  the  receiver  walls  and  the  tem- 
perature potential  between  the  walls  and  the  gas.   Considering 
first  the  denominator,  it  can  be  seen  that  when  the  tempera- 
ture of  the  gas  and  that  of  the  walls  are  the  same,  this  term 
is  zero.   Of  course  in  this  situation,  the  heat  transfer  (the 
numerator  of  the  expression)  is  also  zero  and  thus  as  would 
be  expected  the  convective  heat  transfer  coefficient  is  un- 
defined.  Now,  examining  the  case  where  there  is  a  finite  but 
small  difference  between  the  wall  and  gas  temperatures,  as 
might  be  the  case  in  a  low  mass  flow  rate  charging  process, 
it  is  seen  that  the  equation  for  h  becomes  very  sensitive  to 
errors  in  the  temperature  of  the  gas  or  of  the  receiver  walls. 
It  has  already  been  observed  that  the  expression  for  the  tem- 
perature of  the  walls  is  insensitive  to  small  errors  in  the 
values  of  the  system  variables,  the  gas  pressure  and  mass. 
Unfortunately  this  is  not  the  case  for  the  perfect  gas  rela- 
tion used  to  express  the  gas  temperature.   Even  though  an 
error  in  T  may  seem  negligible  when  considered  on  a  percentaae 
basis,  it  is  not  when  considering  small  differences  between  T 

and  T  .   For  example,  an  error  of  1°R  in  the  temperature  of 

w 

the  gas  will  cut  the  value  for  the  denominator  of  h  in  half 
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if  the  temperature  difference  is  only  2°R.   A  temperature  error 
of  this  magnitude  is  definitely  within  the  experimental  accu- 
racy of  these  tests,  therefore  one  must  exercise  caution  when 
applying  this  equation  to  runs  in  which  the  temperature  poten- 
tial between  the  receiver  walls  and  the  gas  is  small. 

The  numerator  of  equation  (9)  can  be  analyzed  in  much  the 
same  manner.   As  was  the  case  with  the  denominator,  it  can  be 
shown  (by  examining  the  normal  magnitude  of  the  various  terms) 
that  except  for  the  case  of  small  temperature  differences  this 
expression  is  capable  of  absorbing  small  errors  in  P  and  M 
without  producing  a  significant  error  in  h. 

Combining  these  results,  it  is  concluded  that  equation  (9) 
is  insensitive  to  small  error  in  P  and  M  when  the  temperature 
difference  between  the  gas  and  the  tank  walls  is  at  least  of 
the  order  of  10°R  or  greater.   Below  this  value  the  calculated 
values  of  the  heat  transfer  coefficient  may  be  subject  to 
substantial  error. 

The  effect  this  conclusion  has  on  the  confidence  in  the 
results  of  these  tests  is  discussed  in  the  following  sections. 
The  actual  value  of  h  is  used  only  in  the  equation  for  the 
value  of  P*  in  the  isothermal  sink  solution.   The  runs  that 
correspond  to  this  model  have  substantial  temperature  poten- 
tials throughout,  therefore  the  values  of  P*  can  be  calculated 
with  confidence.   In  the  other  two  cases  considered,  the  val- 
ues of  P*  are  not  a  function  of  the  convective  heat  transfer 
coefficient,  but  rather  the  value  of  h  is  used  in  determining 
the  range  of  NTU  in  order  to  determine  the  oartic  alar  solution 
to  be  employed.   The  values  of  NTU  for  these  runs  are  dictated 
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to  a  greater  extent  by  the  mass  flow  rate  than  by  the  convec- 
tive  heat  transfer  coefficient.   For  this  reason,  although  the 
values  of  NTU  may  not  be  exact,  they  can  certainly  be  used  to 
evaluate  the  relationship  between  NTU  and  the  deviation  of  the 
theoretical  P*  values  from  those  actually  observed  in  either 
of  these  cases. 

The  situation  is  less  clear  cut  when  it  comes  to  inter- 
preting the  convective  heat  transfer  coefficient  data.   The 
sensitivity  of  the  expression  for  h  at  low  values  of  the  tem- 
perature potential  places  the  numerical  values  for  some  runs 
in  doubt.   The  low  mass  flow  rate  runs  in  particular  are 
greatly  affected  by  this  uncertainty.   For  these  runs  the  tem- 
perature potential  was  only  on  the  order  of  10°R  for  the  ini- 
tial runs  and  decreased  as  the  charging  pressure  was  progres- 
sively lowered  for  each  successive  set  of  runs.   This  does  not 
mean  that  this  data  is  not  meaningful.   There  is  no  reason  to 
believe  that  the  trends  indicated  by  such  data  should  not  be 
correct;  however,  it  must  be  realized  that  the  numerical  re- 
sults are  subject  to  a  large  range  of  experimental  uncertainty 

B.   COMPARISON  OF  REYNOLDS'  MODELS  WITH  THE  EXPERIMENTAL 

RESULTS 

1.   High  Mass  Flow  Rate  Runs 

As  was  mentioned  earlier,  all  runs  of  this  nature  were 

meant  to  correspond  to  Reynolds'  adiabatic  solution  criteria. 

The  value  for  C*  was  657  in  all  cases  and  the  calculated  val- 
o 

ues  of  NTU  were  well  below  .25  in  each  run.   Reynolds  would 
then  recommend  that  the  solution  to  this  tvoe  of  charcring  be 
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approximated  by  the  assumption  that  there  was  negligible  heat 
transfer  to  the  receiver  walls.   When  a  comparison  of  the  val- 
ues of  the  dimensionless  pressure  as  predicted  by  this  solution 
and  those  observed  experimentally  was  made  the  deviation  be- 
tween the  two  values  was  of  the  order  of  7.5%  of  the  actual 
value  (see  Table  2) ,  with  the  theoretical  value  being  higher 
than  the  actual. 

The  discrepancy  between  these  two  figures  may  be  at  least 
partially  accounted  for  by  noting  that  the  assumption  that 
there  is  no  heat  transfer  is  not  exact.   Since  the  effect  of 
heat  transfer  is  to  lower  the  change  in  internal  energy  of  the 
gas  and  thus  diminish  the  increase  in  the  gas  temperature,  a 
lower  gas  temperature  results  with  a  subsequent  reduction  in 
pressure.   With  this  in  mind  one  expects  the  theoretical  P* 
value  to  be  higher  than  that  actually  observed,  as  is  the  case 
for  this  data. 

The  parameter  NTU  is  a  measure  of  the  thermal  conductance 
(h)  relative  to  the  rate  of  gas  flow  (w~) .   The  value  of  NTU 
should  then  be  related  to  the  exactness  of  the  assumption  of 
no  heat  transfer.   A  connection  between  NTU  and  the  deviation 
of  the  actual  and  theoretical  P*  values  can  be  noted  in  a 
single  run.   NTU  is  directly  proportional  to  h  for  anv  partic- 
ular run,  for  w0  is  a  constant.   An  increase  in  h  for  runs  of 
this  sort  is  accompanied  by  an  increase  in  the  deviation  (Fig. 
5) .   Although  the  quantitative  relationship  between  h  and  wQ 
is  not  clear  from  the  data  for  different  runs,  the  increase  in 
h  and  wQ  is  such  that  their  ratio,  and  hence  NTU,  increases 
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with  increasing  flow.   The  expected  increase  in  the  departure 
of  theory  from  experiment  with  increasing  NTU  is  observed  in 
the  comparisons  (Fig.  8) . 

Reynolds  proposed,  on  theoretical  grounds,  that  for  NTU 
less  than  .25  the  deviation  between  actual  values  of  P*  and 
T*  and  those  predicted  by  the  adiabatic  solution  would  be  less 
than  5%  for  values  of  M*  as  high  as  9 .   In  the  experimental 
data,  this  figure  seemed  valid  at  the  lowest  values  of  NTU 
(-.08)  recorded  and  appeared  slightly  optimistic  at  higher 
values  (10%  deviation  at  NTU  -.12).   In  any  case,  the  accu- 
racy of  these  tests  was  insufficient  to  discriminate  a  5% 
deviation  and  the  criteria  developed  by  Reynolds  is  suffi- 
cient for  engineering  calculations. 

In  conclusion,  it  has  been  seen  that  the  deviation  be- 
tween the  actual  thermodynamic  state  of  a  gas  and  that  pre- 
dicted by  the  adiabatic  solution  for  charging  is  a  function 
of  the  actual  magnitude  of  the  heat  transfer  in  the  system. 
The  actual  temperature  and  pressure  were  always  found  to  be 
lower  than  that  predicted  by  the  adiabatic  solution  for  the 
receiver.   The  parameter  NTU  proposed  by  Reynolds  seems  to  be 
a  good  measure  of  whether  or  not  an  adiabatic  solution  should 
be  attempted.   His  criteria  seems  satisfactory  for  all  but 
the  most  exacting  engineering  calculations  as  well.   There  is 
also  evidence  to  the  fact  that  the  values  of  the  convective 
heat  transfer  coefficient  and  the  mass  flow  rate  are  not  inde- 
pendent in  charging  processes  of  this  nature. 


45 


2.   Low  Mass  Flow  Rate  Runs 

As  was  discussed  in  an  earlier  section,  the  runs  of  this 
type  were  performed  in  order  to  be  used  in  evaluating  Reynolds1 
isothermal  model.   The  value  for  C*  in  all  cases  was  65  7  and 
the  values  of  NTU  ranged  from  1.8  to  14.6.   For  values  of  NTU 
greater  than  7 ,  Reynolds  recommends  the  use  of  the  isothermal 
closed  form  solution  for  charging  processes.   When  the  values 
of  the  dimensionless  pressure,  P*,  as  predicted  by  the  isother- 
mal solution  and  those  observed  experimentally  were  compared 
the  deviation  between  the  two  values  was  somewhere  between  .5% 
and  4.0%  of  the  actual  value  (see  Table  3).   In  all  cases  this 
deviation  proved  to  be  negative;  that  is,  the  actual  pressure 
and  temperature  were  higher  than  those  predicted  in  an  iso- 
thermal solution  to  the  charging  process. 

The  deviation  in  the  isothermal  values  of  P*  and  the 
actual  values  for  these  runs  were  small  and  can  be  accounted 
for  by  investigating  the  assumptions  that  lead  to  the  iso- 
thermal charging  model.   In  this  model  the  thermal  capac- 
itance of  the  receiver  was  assumed  infinite  while  the  internal 
heat  transfer  resistance  was  taken  to  be  zero.   Therefore,  any 
tendency  for  the  incoming  gas  to  heat  up  was  countered  by  an 
immediate  transfer  of  heat  to  the  receiver  walls.   This  con- 
dition was  approached  by  charging  so  slowly  that  the  gas  and 
the  walls,  initially  at  the  same  temperature,  remained  at  that 
state.   This  model  can  not  be  exact  since  the  interfaces  be- 
tween the  gas  and  the  walls  have  a  finite  heat  transfer  resis- 
tance, thus  the  resulting  discrepancies  between  the  actual  and 
theoretical  values  for  P*. 
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It  seems  reasonable  to  predict  in  view  of  the  preceding 
discussion  that  an  increase  in  the  heat  transfer  resistance 
(a  decrease  in  h)  would  lead  to  an  increased  deviation  in  the 
isothermal  and  actual  values  for  P*.   Although  the  values  for 
h  were  subject  to  error  in  these  runs  (see  section  VII-A) ,  a 
definite  downward  trend  throughout  a  single  run  was  noted 
(Fig.  15,  16,  17  and  18).   As  predicted,  this  decrease  in  h 
was  accompanied  by  an  increase  in  the  deviation  of  the  P* 
values  (Fig.  6) . 

An  examination  of  the  effect  of  NTU  on  the  deviation  of 
the  isothermal  and  actual  thermodynamic  state  of  the  gas  was 
then  made.   It  was  noted  that  the  values  of  h  for  different 
runs  remained  constant  within  experimental  scatter  regardless 
of  the  value  of  wQ.   Thus,  unlike  the  adiabatic  case  where 
the  effects  of  changing  wQ  seemed  to  be  offset  by  correspond- 
ing changes  in  h,  the  value  of  NTU  for  this  type  of  charging 
was  inversely  proportional  to  the  value  of  wn.   Once  again 
NTU  proved  to  be  a  good  measure  of  the  deviation,  for  as  w~ 
decreased  the  percentage  deviation  decreased  (Fig.  9) . 

The  criteria  proposed  by  Reynolds  suggests  that  for  values 
of  NTU  larger  than  7  an  isothermal  solution  for  a  charging 
process  will  give  less  than  5%  deviation  from  the  actual  ther- 
modynamic state  of  the  gas.   The  data  from  these  tests  indi- 
cates that  for  values  of  NTU  as  low  as  2  the  deviation  is 
still  less  than  5%,  therefore  it  seems  that  Reynolds'  criteria 
is  slightly  conservative.   Of  course,  it  must  be  pointed  out 
that  the  value  of  C*  for  these  runs  was  65  7  and  this  is  well 
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above  the  value  of  C*  >  40  listed  alonq  with  the  NTU  >  7  cri- 
teria.  Nonetheless,  it  can  safely  be  asserted  that  for  values 
of  C*  above  40  and  values  of  NTU  of  7  or  greater,  the  isother- 
mal solution  to  a  charging  process  will  give  good  accuracy. 

In  conclusion,  it  was  observed  that  the  actual  temperature 
and  pressure  in  a  receiver  being  charged  are  consistently 
higher  than  those  predicted  by  an  isothermal  solution.   The 
parameter  proposed  by  Reynolds  once  again  anpeared  to  be  a  good 
indicator  as  to  whether  or  not  an  assumption  of  isothermal 
charging  can  be  applied  to  a  particular  process.   Reynolds1 
criteria  also  proved  to  be  very  satisfactory  for  engineering 
purposes.   Finally,  evidence  was  observed  that  the  value  of 
the  convective  heat  transfer  coefficient  and  the  mass  flow 
rate  can  be  considered  independent  in  a  charging  process  of 
this  nature. 

3.   Isothermal  Sink  Runs 

The  last  group  of  runs  to  be  evaluated  were  those  made 
with  the  test  receiver  immersed  in  an  ice  bath.   These  runs 
were  performed  at  intermediate  values  of  NTU  so  as  to  corre- 
spond to  Reynolds'  isothermal  sink  closed  form  solution.   The 
only  criterion  specified  for  the  use  of  this  method  is  that 
the  dimensionless  capacitance  term  be  greater  than  40.   In  all 
the  runs  of  this  nature  the  value  for  C£  was  626,  therefore 
this  criteria  was  easily  satisfied.   Of  course,  Reynolds  im- 
plies that  the  value  of  NTU  dictating  an  isothermal  sink  solu- 
tion (to  approximate  the  actual  conditions)  be  between  .25  and 
7,  for  if  not  the  regions  encountered  prescribe  either  the 
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adiabatic  or  the  isothermal  solution.   When  the  values  of  the 
dimensionless  pressure  as  predicted  by  the  isothermal  sink 
solution  and  those  observed  experimentally  were  compared,  the 
deviation  between  the  two  values  was  somewhere  between  8%  and 
14.5%  of  the  experimental  value. 

It  should  be  noted  here  that  the  isothermal  sink  solution 
is  an  intermediate  case  between  the  two  extremes  of  an  adia- 
batic solution  and  the  isothermal  solution.   Therefore,  the 
thermodynamic  state  of  the  gas  must  be  expressed  as  a  func- 
tion of  the  heat  transfer  coefficient  and  the  mass  flow  rate 
in  the  closed  form  solution  (see  equations  (4)  and  (5)). 

As  before,  physical  interpretations  were  souqht  for  the 
deviation  between  the  isothermal  sink  solution  for  the  thermo- 
dynamic state  of  the  gas  and  that  actually  observed.   The 
assumption  on  which  this  model  was  built  was  that  the  temper- 
ature of  the  thermal  capacitance  remained  a  constant  even 
after  absorbing  the  energy  released  by  the  gas  in  the  receiver, 
If  the  comparatively  small  amounts  of  heat  transferred  to  the 
receiver  walls  were  uniformly  distributed  throughout  the  very 
large  thermal  capacitance,  this  assumption  might  be  quite 
accurate.   However,  as  was  seen  by  an  investigation  using  the 
heat  conduction  equation  for  a  semi-infinite  solid  (see  sec- 
tion V-D)  most  of  the  internal  energv  (as  indicated  by  the 
temperature  distribution)  never  gets  any  farther  than  the 
first  few  tenths  of  an  inch  in  a  charaing  process  of  this  na- 
ture.  For  this  reason  the  temperature  at  the  wall  surface 
does  not  remain  constant  but  rather  increases  slightly  even 
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though  the  effective  thermal  capacitance  of  the  receiver  is 
still  extremely  large.   Bearing  this  in  mind,  the  actual  tem- 
perature potential  between  the  inside  surface  of  the  receiver 
walls  and  the  gas  is  overestimated  by  the  isothermal  sink 
solution.   Therefore,  the  amount  of  heat  transferred  out  of 
the  gas  and  the  accompanying  temperature  and  pressure  drops 
are  slightly  exaggerated  in  this  closed  form  solution  (Fig. 
7). 

To  begin  an  evaluation  of  the  parameters  that  affect  the 
deviation  between  theory  and  experiment,  single  runs,  in 
which  the  mass  flow  rates  were  constant,  were  examined.   For 
the  low  mass  flow  runs  there  seemed  to  be  only  a  very  gradual 
increase  in  the  value  of  h  as  the  run  proceeds.   At  higher 
values  of  the  flow  rate  the  increases  in  h  during  a  run  were 
more  substantial  as  were  their  magnitudes  (Figs.  19-21). 
Using  both  the  low  mass  flow  runs  (almost  constant  h)  and  the 
high  flow  rate  runs  (increasing  h)  it  was  possible  to  examine 
the  effect  of  h  on  the  deviation  between  an  actual  run  and 
the  thermodynamic  state  predicted  for  this  run  by  the  iso- 
thermal sink  solution.   In  analyzing  the  data  (Figs.  7  and 
10)  the  conclusion  was  reached  that  the  deviation  was  appar- 
ently not  closely  dependent  on  the  value  of  h.   In  both  the 
low  flow  rate  runs  and  the  high  flow  rate  runs  the  deviation 
followed  the  same  pattern  throughout  a  run  and  likewise  there 
was  no  discernible  connection  between  the  magnitude  of  the 
deviation  and  the  values  for  h.   In  any  case  it  is  probable 
that  the  effect,  if  any,  of  the  value  of  h  on  the  deviation 
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of  experimental  and  isothermal  sink  values  is  masked  by  the 
use  of  NTU  =  hA/w0c   as  a  parameter  in  equations  (4)  and  (5) . 

For  the  two  extreme  solutions  (adiabatic  and  isothermal) 
presented  in  the  previous  sections ,  their  accuracy  increased 
as  the  appropriate  extreme  condition  was  approached.   This 
might  lead  one  to  believe  that  the  intermediate  ranges  of  NTU 
for  the  isothermal  sink  solution  now  being  considered  would 
most  closely  correlate  to  the  actual  thermodynamic  state  of 
the  gas  during  charging.   Since  the  data  does  not  cover  the 
entire  range  of  NTU  between  .25  and  7,  no  conclusive  statement 
can  be  made  in  this  regard.   Figure  10,  however,  supports  the 
notion  that  as  NTU  approaches  an  intermediate  value  of  say  3.5 
the  deviation  decreases.   One  thing  can  be  asserted,  and  that 
is  that,  due  to  the  connection  between  wQ  and  h  in  this  type 
of  charging,  the  values  of  NTU  do  not  vary  greatly.   Thus, 
regardless  of  the  connection  between  NTU  and  the  deviation, 
it  seems  certain  that  the  deviation  between  the  isothermal 
sink  solution  and  that  actually  observed  is  approximately 
constant  for  a  particular  system  as  long  as  the  flow  rates 
are  of  the  same  order  of  magnitude. 

In  conclusion,  it  should  be  noted  that  the  isothermal 
sink  closed  form  solution  is  a  function  of  the  convective 
heat  transfer  coefficient  and  the  mass  flow  rate,  and  thus 
differs  from  the  previous  two  closed  form  solutions.   The 
deviation  between  the  thermodynamic  state  of  the  gas  as  pre- 
dicted by  the  analytical  solution  and  that  observed  experi- 
mentally was  between  8%  and  15%  of  the  experimental  value. 
The  closed  form  solution  values  were  consistently  lower  than 
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those  found  ±n  the  tests.   This  discrepancy  can  be  explained 
by  the  fact  that  the  predicted  value  of  the  temperature  of 
the  thermal  capacitance  in  Reynolds '   solution  is  lower  than 
that  actually  existing.   The  heat  transfer  and  temperature 
drop  are  subsequently  overestimated  in  this  solution.   For 
runs  of  this  nature  it  was  also  observed  that  the  values  of 
h  increase  throughout  a  run  and  that  not  only  the  maanitude 
of  this  increase  but  the  values  of  h  are  increased  when  the 
mass  flow  rate  is  increased.   There  seemed  to  be  little  or  no 
dependence  of  the  accuracy  of  the  isothermal  sink  assumption 
on  the  values  for  h  when  wQ  is  a  constant.   No  conclusive 
statement  can  be  made  as  to  the  connection  between  NTU  and 
the  deviation  of  the  isothermal  sink  solution  and  the  actual 
values ,  although  there  is  evidence  to  support  a  prediction 
that  this  solution  is  most  accurate  at  intermediate  values  of 
NTU.   Due  to  the  apparent  connection  between  w~  and  h  the 
values  of  NTU  do  not  change  substantially  for  runs  in  this 
region  nor  do  the  percentage  deviation  in  the  values  for  the 
thermodynamic  state  of  the  gas.   Lastly,  for  this  particular 
experiment  the  discrepancy  between  the  solution  as  predicted 
by  an  isothermal  sink  solution  and  that  observed  experimen- 
tally, was  slightly  higher  than  that  predicted  by  Reynolds 
even  though  the  value  of  C*  was  well  above  the  criteria  value 
of  40. 
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C.   DISCUSSION  OF  THE  EXPERIMENTALLY  OBTAINED  HEAT  TRANSFER 

COEFFICIENTS 

Having  empirically  obtained  values  for  the  convective  heat 
transfer  coefficient  throughout  the  various  runs,  an  attempt 
was  made  to  obtain  at  least  a  qualitative  explanation  of  the 
behavior  of  these  values.   The  charging  situation  was  sub- 
stantially different  from  those  usually  described  in  convection 
theories.   The  pressure  of  the  system  is  usually  not  varied 
in  such  theories,  nor  is  the  direction  and  magnitude  of  the 
flow  velocity  a  function  of  variables  other  than  position. 
Since  any  attempt  to  arrive  at  a  theoretical  expression  for 
the  convective  heat  transfer  coefficient  in  a  receiver  being 
charged  would  be  completely  beyond  the  scope  of  this  paper, 
no  attempt  was  made.   Some  insight  was  gained,  however,  in 
considering  the  possible  non-dimensional  parameters  that 
might  affect  an  expression  for  the  convective  heat  transfer 
in  a  charging  situation. 

From  the  general  equations  of  motion  and  enercry ,  an 
average  convective  heat  transfer  coefficient  evaluated  over 
the  total  surface  area  of  a  body  (this  eliminates  any  spatial 
dependence)  can  be  expressed  as  a  function  of  the  dimension- 
less  parameters  as  listed  below: 

NU  =  hL/k  =  f(PR,  ER,  GR ,  RE) 
An  investigation  of  these  parameters  was  th^n  made  in  order 
to  come  to  a  better  understanding  of  the  causes  for  the  behav- 
ior of  h  in  the  experimental  data.   First  to  be  considered 
was  the  Prandtl  number.   Defined  as 

PR  =  yc  /k 
P 
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this  dimensionless  number  can  be  considered  constant  and  of 
the  order  of  1  for  solutions  involvinq  perfect  gases.   Having 
concluded  that  the  effects  of  the  Prandtl  number  need  not  be 
considered,  the  Eckert  number  was  then  considered.   This  non- 
dimensional  parameter  expresses  the  ratio  of  the  inertia  terms 
to  the  enthalpy  terms  and  is  defined  as 

ER  =  U2/c  AT 
P 

where  U  is  the  characteristic  velocity  of  the  flow  and  c  AT 

p 

is  the  change  in  enthalpy  across  the  boundary  layer.   In  the 
charging  processes  carried  out  for  this  paper  the  average  gas 
velocities  were  small  as  compared  to  the  enthalpy,  therefore 
the  effect  of  the  Eckert  number  is  considered  negligible. 

As  has  been  pointed  out,  a  common  practice  is  to  take  the 
Grashof  number  as  the  governing  parameter  for  the  heat  trans- 
fer processes  of  the  kind  under  consideration  here.   Defined 
as  , 

GR  =  gL3SATp2/U2 

this  term  is  the  ratio  of  the  buoyancy  to  the  viscous  effects. 
The  buoyancy  of  the  fluid  aids  the  heat  transfer  by  supporting 
the  motion  of  the  fluid  at  the  surface  of  the  receiver  walls. 
Noting  this,  the  Grashof  number  must  be  considered  significant 
in  determining  the  convective  heat  transfer  coefficient  in  a 
charging  process  of  the  sort  presented  here.   Since  consider- 
able motion  is  also  present  in  the  tank  due  to  the  stirring 
action  of  the  injection,  the  Reynolds  number,  defined  as 

re  =  up£/y 
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must  also  be  considered.   In  the  above  expression,  U  is  again 
a  characteristic  velocity  of  the  flow  and  I    a  characteristic 
length  for  the  flow.   This  dimensionless  number  can  be  seen 
to  be  a  measure  of  the  relative  importance  of  the  inertia  and 
viscous  effects  and  often  occurs  in  the  analysis  of  problems 
in  forced  convection. 

Investigating  the  variables  comprising  these  non-dimen- 
sional parameters  leads  to  some  understanding  as  to  their 
magnitude  and  behavior  for  charging  processes.   For  small 
temperature  changes  the  density  can  be  considered  proportional 
to  the  pressure.   If  the  temperature  changes  are  significant 
the  density  is  both  proportional  to  the  pressure  and  inversely 
proportional  to  the  temperature.   The  viscosity  is  almost 
independent  of  pressure  and  increases  slightly  with  tempera- 
ture.  For  situations  of  an  isothermal  nature,  it  is  then  ex- 
pected that  the  variables  in  the  Reynolds  number  are  the 
velocity,  the  density  (which  is  proportional  to  the  pressure) 
and  perhaps,  the  characteristic  length  of  the  fluid  flow.   If 
the  temperature  changes  throughout  a  process  are  large,  once 
again  the  Reynolds  number  is  proportional  to  the  velocity,  the 
density  (now  proportional  to  the  pressure  and  inversely  pro- 
portional to  the  temperature),  the  characteristic  length,  and 
also  inversely  proportional  to  the  viscosity  (which  is  propor- 
tional to  the  temperature) .   Now  applying  the  same  sort  of 
reasoning  to  the  Grashof  number,  the  effect  of  y  as  a  variable 
in  a  moderate  temperature  change  system  can  acrain  be  eliminated 
along  with  the  gravitational  constant,  g.   The  density  behaves 
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in  the  manner  described  in  the  Reynolds  number  discussion. 
The  coefficient  of  volumetric  exDansion,  $,  is  defined  as 

3  =  -l(iP) 
P     P  3T;p 

and  varies  as  the  inverse  of  the  temperature  for  a  perfect 
gas.   The  term  AT  in  the  Grashof  number  is  the  temperature 
change  across  the  boundary  layer  and  L  is  the  characteristic 
length  of  the  system,  a  constant. 

In  preparation  for  applying  what  was  learned  about  the 
factors  affecting  the  heat  transfer  coefficient  to  the  behav- 
ior observed  in  the  experimental  runs,  the  velocity  variations 
in  the  test  receiver  were  examined.   When  the  inlet  valve  for 
the  receiver  was  first  opened  a  pressure  ratio  (the  receiver 
back  pressure/  the  stagnation  pressure  upstream  of  the  flow 
meter)  considerably  smaller  than  the  critical  oressure  ratio 
existed  across  the  flow  metering  device.   Sonic  flow  was 
established  at  the  throat  of  the  inlet  to  the  tank.   Therefore, 
fluid  entering  at  a  pressure  corresponding  t_o  the  sonic  condi- 
tions was  forced  to  adjust  to  the  lower  back  pressure  in  the 
tank.   Regardless  of  the  form  of  this  adjustment  it  must  have 
been  coupled  to  a  substantial  increase  in  the  local  velocity. 
As  the  back  pressure  in  the  receiver  increased,  the  pressure 
adjustment  became  less  severe  and  thus  the  local  velocity  de- 
creased throughout  a  charging  process.   Here  the  term  "local" 
applies  to  the  conditions  in  the  tank  near  the  point  of 
injection. 

For  the  low  mass  flow  rate  runs  the  temperatu  e  potential 
between  the  gas  aad  one  wails,  AT,  was  very  small.   For  this 
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reason  it  was  felt  that  the  Grashof  number  was  not  as  signif- 
icant as  the  Reynolds  number  throughout  these  runs.   Having 
seen  that  during  a  charge  the  velocity  near  the  injector  de- 
creased, it  was  thought  that  perhaos  the  characteristic  length 
over  which  the  incoming  jet  acts  was  also  decreased.   There- 
fore, although  the  density  increased  due  to  the  increase  in 
pressure,  it  seemed  reasonable  to  assume  that  the  Reynolds 
number  decreased  slightly  throughout  a  run.   A  droo  in  the 
Nusselt  number,  and  thus  the  convective  heat  transfer  coeffi- 
cient, was  predicted.   In  other  words,  forced  convection  gives 
way  to  free  convection  as  the  run  slowly  proceeds.   Figures 
15,  16,  17  and  18  seem  to  support  a  conclusion  of  this  sort. 
This  notion  also  is  supported  in  a  paper  by  Ulrich  et.  al. 
[Ref.  6]  in  which  the  mode  of  heat  transfer  following  gas 
injection  was  observed  to  proceed  from  forced  to  free  convec- 
tion.  For  the  high  mass  flow  rate  and  isothermal  sink  runs, 
the  temperature  potential  was  much  greater  than  that  for  the 
low  mass  flow  rate  runs,  thus  the  Grashof  number  was  exoected 
to  be  more  significant  in  these  processes.   Since  the  Grashof 
number  is  greatly  affected  by  a  pressure  increase  as  dictated 
by  the  p2  term,  it  was  felt  that  the  Nusselt  number  would  in- 
crease throughout  a  run.   Once  again  these  assumptions  seem 
to  be  supported  by  the  data,  see  figures  11,  12,  13  and  14  as 
well  as  19,  20  and  21. 

As  to  the  actual  magnitude  of  the  convective  heat  transfer 
coefficient  as  indicated  by  NTU,  the  data  was  hard  to  corre- 
late.  As  seemed  consistent  with  the  conclusions  drawn  above, 
the  values  of  the  Nusselt  and  Grashof  numbers  for  the  low  mass 
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flow  rate  runs  were  lower  than  those  for  the  high  mass  flow 
rate  and  isothermal  sink  cases.   However,  with  the  exception 
of  the  set  of  runs  at  a  very  low  mass  flow  rate  (perhaps  cor- 
responding to  the  low  AT  discussion  above) ,  the  isothermal 
sink  values  were  above  those  for  the  high  mass  flow  rate  runs. 
One  possible  explanation  for  this  was  the  fact  that  the  iso- 
thermal sink  runs  were  recorded  in  1  second  intervals  for  run 
times  of  7  to  9  seconds.   The  high  mass  flow  rate  runs,  on 
the  other  hand,  were  only  of  1  second  duration  with  the  points 
measured  at  1/10  second  intervals.   Therefore  it  might  be  that 
for  runs  of  larger  duration,  the  mixing  motion  caused  by  the 
impinging  jet  could  spread  more  thoroughly  throughout  the  re- 
ceiver.  These  types  of  forced  convection  phenomena  could  not 
be  adequately  described  on  the  basis  of  the  results  of  these 
tests.   In  any  case,  the  values  calculated  in  these  tests  are 
at  most  1  order  of  magnitude  above  those  predicted  by  free 
convection  theories  (Figs.  11-21).   Therefore  if  no  better 
estimate   of  the  heat  transfer  coefficient  can  be  obtained, 
a  value  based  on  these  theories  could  be  used  for  processes 
of  the  nature  discussed  here. 

The  preceding  discussion  was  based  mainly  on  physical 
intuition  and  is  only  crudely  substantiated  by  the  few  experi- 
mental runs  performed  in  this  undertaking.   As  will  be 
repeated  in  the  recommendations,  a  much  more  extensive  and 
controlled  study  is  needed  in  this  area  to  draw  any  conclu- 
sive relationships  between  the  convective  heat  transfer  co- 
efficient and  the  many  variables  in  a  system  being  charged. 
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VIII.   CONCLUSIONS 

The  general  conclusions  to  be  drawn  from  this  investiga- 
tion are  summarized  as  follows. 

(a)  Excellent  quantitative  agreement  between  the 
experimental  values  and  those  predicted  by  Reynolds 
[Ref.  1]  has  been  obtained  for  the  thermodvnamic 
state  of  a  gas  in  a  receiver  being  charged. 

(b)  The  parameter  NTU  prescribed  by  Revnolds  is  a  good 
indicator  for  describing  the  extent  of  heat  transfer 
in  a  receiver  being  charged. 

(c)  The  technique  developed  in  this  paper  for  experi- 
mentally determining  the  convective  heat  transfer 
coefficient  may  be  confidently  used  in  conjunction 
with  Reynolds'  closed  form  solutions  if  the  tem- 
perature potential  between  the  receiver  walls  and 
the  charging  gas  is  of  an  order  of  10°R  or  larger. 

(d)  The  convective  heat  transfer  coefficient  is  related 
to  the  initial  mass  flow  rate  of  the  entering  gas 
for  all  but  the  extremely  slow  charcring  rates 

(wQ  <  .003  lbm/sec). 

(e)  The  heat  transfer  coefficient  for  charging  of  the 
nature  described  in  this  paper  can  be  estimated 
within  one  order  of  magnitude  by  free  convection 
theories.   The  heat  transfer  coefficient  thus 
obtained  will  be  equal  to  or  less  than  the  actual 
value.   The  effects  of  errors  in  this  estimate 
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will  be  small  provided  that  the  NTU  is  sufficiently 
small  (high  injection  rates). 

IX.   RECOMMENDATIONS 

(1)  It  is  recommended  that  an  extensive  and  controlled  studv 
of  the  effects  of  the  various  system  parameters  on  the 
convective  heat  transfer  coefficient  during  charcring  be 
made  with  emphasis  placed  on  the  development  of  an 
empirical  formula  for  use  in  this  area. 

(2)  A  possible  extension  of  the  technique  for  experimentally 
determining  the  convective  heat  transfer  coefficient 
developed  in  this  paper  could  include  the  use  of  analog 
computer  methods  to  directly  analyze  the  pressure  and 
temperature  signals  from  probes  in  the  vessel. 

(3)  In  any  further  investigation  of  the  type  discussed  in 
this  paper,  the  following  suggestions  are  made.   A  more 
exact  method  of  measuring  the  mass  flow  rate  should  be 
developed.   Care  should  also  be  taken  to  eliminate  any 
possible  transient  effects  in  establishing  flows  used  in 
experiments  of  this  nature.   It  is  also  recommended  that 
hot  gases  be  used  in  order  to  obtain  substantial  temper- 
ature potentials  between  the  receiver  walls  and  the  qas. 
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APPENDIX  A 

Nomenclature  of  Reference  1 

English  Letter  Symbols 

A        Area,  ft2 

C        Thermal  Capacitance  of  receiver  shell,  Btu/°R 

c        Specific  heat  at  constant  pressure,  Btu/(lb°R) 

c        Specific  heat  at  constant  volume,  Btu/(lb°R) 

h  Specific  enthalpy,  Btu/lb 

h        Unit  heat  transfer  convective  conductance, 
Btu/(hr  ft2oR) 

k  Thermal  conductivity,  Btu/(hr  ft2oR/ft) 

M  Mass,  lb 

P  Pressure,  lb/ft2 

q  Heat  transfer  rate,  Btu/hr 

T  Absolute  temperature,  °R 

U  Total  internal  energy,  Btu 

u  Specific  internal  energy,  Btu/lb 

V  Volume,  ft3 

w  Mass  flow  rate ,  lb/hr 

Greek  Letter  Symbols 

0  Time,  hr 

y  Viscosity,  lb/(hr  ft) 

p  Density,  lb/ft3 

A  Denotes  a  difference 
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Nondimensional   Grouping 


c* 

C    /(M„c    ) 
&       0    v 

k 

c   /c 
p       V 

M* 

M/MQ 

NTU  l/(RiCvw0)    =     (hA)i/cvwQ 

1/(R   c   wn)    =    (hA)    /c  wn 

00  V  0  oo   V  0 


NTU 

oo 

V(RW 

p* 

p/p0 

T* 

T/TQ 

T* 

T1/TQ 

T* 
c 

Tc/T0 

T* 

00 

Too/T0 

W* 

w/wQ 

Subscrd 

.pts 

0  Refers  to  initial  conditions 

1  Refers  to  inlet  state 
c  Refers  to  capacitance 

i  Refers  to  inside  receiver 

00  Refers  to  environmental  conditions  outside  receiver 
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APPENDIX  B 

Derivation  of  General  Differential  Equation  and  Closed  Form 

Solutions  for  Charging 
Assumptions  and  Idealizations 

(1)  The  walls  of  the  receiver  are  at  a  uniform  temperature 
throughout. 

(2)  The  heat  transfer  resistances  both  inside  and  outside 
are  constant  and  uniform  throughout  their  respective 
surface  areas. 

(3)  Perfect  mixing  of  the  injected  gas  and  the  gas  in  the 
receiver. 

(4)  The  working  fluid  is  a  thermallv  and  calorificallv 
perfect  gas. 

(5)  The  specific  heat  of  the  capacitance  is  a  constant. 

(6)  The  kinetic  energy  of  the  gas  in  the  receiver  may  be 
neglected. 

General  Differential  Equation  for  Charging 


/////////////// 


Z//////A 


dM 


///////// 
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An  energy  balance  on  the  receiver,  using  h,    as  the  stagnation 
enthalpy  of  the  inlet  gas  gives , 

,   dM    dU 

hi   de  =  de  +  qi  • 

An  energy  balance  on  the  receiver  walls  yields , 

dUc 

qi  =  de2  +  q-  • 

The  heat  transfer  rates  can  be  defined  in  terms  of  the  temper- 
ature potential  and  the  heat  transfer  resistance  as  follows: 

T-T  T  -T 

_  c  c   °° 

qi  :    R.   '    q~  ::   R 
1  °° 

The  change  in  the  total  internal  energy  of  the  gas  in  the 

receiver  may  be  written  in  terms  of  the  mass  and  the  specific 

internal  energy  as , 

dU   ..du     dM 

de  =  M3e  +  ude  ■ 

The  rate  of  change  of  energy  in  the  capacitance  mav  be  expres- 
sed as  , 

dU    d(C  T  )      dT 
c       c  c   _    __c 

de  ::    de   "  c  de 

Thus,  the  energy  balance  on  the  gas  may  be  written  as, 

T-T 
,   dM    ..du     dM  ,     c 

hi   de  "  Mde  +  ude  +  ~r~ 

i 


or 


and  thus , 


l%        ^   ,.d(cvT>    T"Tc 

(Vu)de  "  M— ae—  +  -r— 


T-T 

dM       dT       c 
(c  T,-c  T  ^  =  c  M^  +  -^-°- 
p  1   v   d6     v  de     R. 
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Note  here  that  T,  is  the  stagnation  temperature  of  the  inlet 
gas.   Therefore,  dividing  by  c   gives, 

(kT  _T)^i  -  MdT    ^£ 
(KT1  T'dT  ~  Md0  +  HT^   ' 

1 

In  a  similar  fashion  the  energy  balance  on  the  capacitance 

can  be  written  as , 

T-T        dT     T  -T 
°  =  C   —2  +   °   °° 


R.      c  d6       R 

l  °° 


then  dividing  by  C  , 


dT           T-T  T    -T 

c         c  c      °° 

d0      :     R.C  '     R  C 

1    c  °°    c 


Therefore , 


dT  _  T  T      t  . 

_£  =      T         +         °°      _    _£(L_  +    1_  ) 
d6  R.C  R   C  C    lR.         R 

1     C  o°     C  CI  °o 

The    instantaneous   mass    flow   rate,    w,    mav  be   defined   as 

dM 
w  =    dT 


and  we  can  write, 


as  well  as, 


dT    dT  dM     dT 
d9  '  dM  d6    WdM 


dT     dT   ,„     dT 

c      c  dM     c 

d6     dM   d6    WdM 
The  equation  for  the  energy  of  the  gas  may  then  be  written  as, 

1  V 


or 


T-T 

wMd¥  +  r-5^  +  (T"kTi>w  =  °  •  (1) 

i  v 
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In  a  similar  manner  the  energy  equation  of  the  capacitance 

can  be  written  as 

dT  T       T   -     . 
c  _   T        °o    -  _c  ,1_    1  v 

WdM   "'  R.C     RC      ClR.R; 

1  C      °°  C        C    1      oo 

or 

/Tc    ^£,1__    1_,  _  _T_  _  J^_  =  0.       (2) 

dM     C   R.    R  '  R.C     R  C 

C   1     °°      1  c     °°  c 

Expressing  equation  (1)  in  dimensionless  form  by  dividing 

through  by  T_  and  noting  that 

-     M*  — £_ 
1  HFF 


M— : 

dM 

M           d 

we  obtain, 

T*  —  T* 

wM*S  +  s +  (T*-kT*)w  =  0  . 

dM*    R.c  1 

l  v 

Now  dividing  by  wQ  and  noting  the  definition  of  NTU,  we  ob- 
tain the  dimensionless  form, 

dT* 
W*M*HM~*  +  (T*-TJ)NTU  +  (T*-kT*)w*  =  0  .      (la) 

In  the  same  way,  equation  (2)  can  be  non-dimensionalized  by 

dividing  through  by  TQw0/MQ  so  that, 

dT*         T*MA    ,  ,  T*MA  T*Mn 

* C  c    0  ,1  1    v  0  °°    0  n 

W    dM*         C   wnlR.  R    '  R.wnC  R   C   wA 

cOi  °°  iOc  °°    c    0 

and  noting  that  NTU/C*  is  equal  to  Mn/R^wnC^  the  following 

is  obtained, 


dT*  NTU 

w*^T7?-  +    T*(NTU+NTU    )     -    S^T*    -    t^t-^T*    =    0    .  (2a) 

dM*  c  °°  Ci  C*         °° 

These  two  dimensionless  equations  can  be  combined  to  vield  a 
single  equation  giving  T*  as  a  function  of  M*  by  differentiation 
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of  equation  (la)  with  respect  to  M*  and  suitable  arithmetic 

manipulations  represented  below: 

dT* 
w*M*3m~*~  +  (T*-T£)NTU  +  (T*-kT*)w*  =  0  .      (la) 

Differentiation  with  respect  to  M*  yields, 

dT* 

rearranging,  the  following  is  obtained, 

dT *        d2T* 

dTJ   w«(2§fr+  m*  ygr)    „, 

dM*  NTU  +  dM* 

dM^(T  +M  ^    W  kTldM*-    . 
NTU  NTU 

Substituting  this  expression  for  dT*/dM*  in  equation  (2a) , 
*2  *d2T*    0  *2dT*  *dw* ,_* , „*dT* 

W*  m*  ■...  ■  m      2W*  -3TTT        Jm*     W*3r75r(T*+M*^-r5-) 

dM*2        dM*     *dT*      dM*       dM* 


NTU         NTU        dM*  NTU 

W*kT*^-    T*  „_.,      NTU 

idM*_    _£(NTU+NTU  )  _  NTU °°  * 

NTU       C*  ^1U™-LU00'    r*        C*   oo 

Multiplying  through  by  NTU  and  combining  gives, 

d2T*  rlw*  r?T* 

w*m*Hm*^  +  W*<2w*+NTU+Mi^>ik 

,    *         .TmTT2  NTU   NTU  ,    * 

+    (w*^-  -   NTU    )T*    -  T*    _   w*kT*dw*_      (3) 

+     [W    dM*         ^^  C*  ioo         W    K1ldM*        <J' 

+    ^(NTU+NTU    )NTU    =    0 
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Now  solving  (la)  for  T*, 

c 


T*  = 
c 


dT* 
w*M*^r-  +  T*NTU+T*w*-kT*w* 
dM* 1 

NTU 


Therefore  the  term, 


T* 
c 


(NTU+NTU  )NTU 


in   the    above   expression   can  be  written   as, 

NTU+NTU^  * 

W*M*H^( C* )     +    ^r^NTU+NTUj 

w*T*                                  kT.w 
+    )~     (NTU+NTU    ) ^—(NTU+NTU    ) 

Substituting  this  in  equation  (3)  we  arrive  at  the  general 
differential  equation  for  charging. 


d2T*  rlw* 

W*M*dW  +  w*[2w*+NTU+M*^  +  ( 


NTU+NTU 

c 

— c* — 


-)  M*  ] 


dT* 
dM*" 


,  *    ,                  NTU  NTU 
+  [w*§^  +  ^(NTU+NTU)  w*  + in-   "IT* 


dM" 


0 


^0 


-  kT*w* 


dw' 


NTU+NTU 


1   dM* 


c* w*kTl  " 


NTU  NTU 

c 

— cr* — 


T*  =  0 


Closed  Form  Solution  for  Adiabatic  Charging 


i  i 

/  / 

//////"A  dU   =    d(Mu)  ^ 

/ 


////////  /J 


M,    T 


/ 
/ 
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Since  there  are  no  heat  transfer  rates  to  consider  in  an 
adiabatic  solution  it  is  not  necessary  to  specify  the  mass 
flow  rate.   An  energy  balance  on  the  gas  in  the  receiver, 
using  stagnation  functions,  gives 

k,  dM  =  d(mu)  =  Mdu  +  udM 
or 


(T^-uJdM  =  Mdu 


Then  one  can  write, 


(c  T..-C  T)  dM  =  Mc  dT 
p  1   v  v 

Dividing  by  c   gives, 

(kTx-T) dM  =  Mdt   . 

In  order  to  get  this  equation  in  dimensionless  form,  we 
divide  through  by  TfiM   thus  obtaining 

(kT*-T*)dM*  =  M*dT* 

Now  separating  variables  and  integrating  this  eauation  from 
the  initial  conditions  (T  =  Tn,  M  =  Mn)  to  a  state  at  a  later 


0 


0' 


time  (T  =  T,  M  =  M) , 


T_ 
T, 


dT* 
(kT*-T*) 


M_ 
M( 

M, 


dM* 

M* 


This  gives , 

-log(k,T*-T*)  +  log(kT*-l)  =  logM5 

Rearranging,  this  equation  reduces  to 


(kT*-l) 

M* 


=  kT*  -  T* 
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Solving  for  T*,  the  final  result  is  obtained, 


T*  =  kT*  - 


(kT*-l) 


Thus  having  derived  a  relation  for  T*  as  a  function  of  M*, 
the  perfect  gas  relation  P*  =  M*T*  is  used  to  obtain, 

P*  =  kT*(M*-l)  +  1   . 

Closed  Form  Solution  for  Isothermal  Charging 


//////////////// 


/ 
/ 

/ 

T    = 

■ 

= .  T  ■ 
c 

■'■}■'■. 

////////  /yy/y/ 


Since  T  =  constant  =  Tn   =   T  ,  then 

0     c 

T 
10 

Therefore  the  isothermal  pressure  mass  relationship  can  be 

obtained  directly  from  the  perfect  gas  law,  and  is 

P*  =  M* 


70 


Closed   Form   Solutions    for   Charging   at   Constant   Mass    Flow  with 
Heat   Transfer   to    an   Isothermal    Sink 


//  /  //  /  /  /  /  / 


dU 

de 


d(Mu) 

de 


X 
X 


x. 


Vo^ldl 


7777777777 


M,    T 


R, 


///// 


X 


zzz 


*c 


An   energy  balance   on   the  gas    using   stagnation    functions   gives, 


,     dM  d(Mu) 

}l^■AK  =  qn-   +  — aK~ 


dM        du.. 
=   q,    +    u-^  +    -^M      . 


By  definition  of  the  heat  transfer  resistance  R. , 

T-T 


qi  "  — 


and  using  w  =  dM/d9 ,  then 

T-T 


c  T,wn  =   _  C  +  R~   +  WnC  T 
p  1  0     R.      d0     0  v 


or 


T-T 


cptiwo-=  — 


c       dT 
-  +  Mc  -r^  +  wnc  T 
vde     0  v 


Using  the  following  relation, 

dT    dT  dM 


dT 


dB    dM  d9    W0dM 
and  dividing  by  wn  and  rearranging  the  exoression  becomes, 

T-T 

Hf  r> 

t-G   T.  +  c  M~  +  c  T  + -  =  0   . 

pi     v  dM     v     R>w» 
r  l  0 
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Now  dividing  by  c   to  obtain 

T-T 
dT  c 

m—   +  T  -  kT.  +  = =  0   . 

dM  1    R.wnc 

l  0  v 

Note  that  all  functions  refer  to  stagnation  conditions.   In 

order  to  get  this  equation  in  dimensionless  form  it  is  divided 

by  TQ  so  that, 

£L  |T  +  T*  _  kT*  +  (T*-T*) 


Since 


Tn  dM  1         c'R.wAc 

0  i  0  v 


M   dT    M      0     A,*dT*    .  „__      1 
7Z—   jT7  =  n n —  -  M*-3T7=r  and  NTU  =  = 

T0  m        M0  d(^->      ^  Riw0Cv 

M0 


this  expression  can  be  written  as, 

M*^t*  +  (T*-kT*)  +  (T*-T*)NTU  =  0   . 
dM*  1  c 

In  order  to  obtain  the  closed  form  solution  the  variables 

are  separated  and  the  expression  integrated  from  the  initial 

conditions  (T  =  TQJ    M  =  MQ)  to  a  state  at  a  later  time  (T  =  T, 


M  =  M) 


T  M 

T0  dT*^ (  M0  dM* 

[T*  (1+NTU)  -kT*-NTUT*]  ==  ~\    M*~ 

To  Mo 


This  gives 


(1+NTU)  log [-kT*-NTUT*+ (1+NTU) T*] 
-  (1+ntu)  !og [~kT*-NTUT*+ (1+NTU) ]  =  -logM' 
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or 

-kT*  -  NTUT*  +  (1+NTU)T*  =  m*~(1+NTU)  [-kT, -NTUT*+ (1+NTU)  ]  . 

Solving  this  expression  for  T*  as  a  function  of  M*  the  follow- 
ing is  obtained: 


T*  = 


kT*+NTUT*- (kT*+NTUT*-l-NTU) M* 
1      C     1      c 


-  (1+NTU) 


(1+NTU) 


Using  the  perfect  gas  law, 


P*  = 


(kT*+NTUT*)M*-(kT*-l-NTU+NTUT*)M* 
1 c        1 c 

1+NTU 


-NTU 


Closed  Form  Solution  for  Charging  at  Constant  Mass  Flow  with 
Inside  Resistance  Negligible 


///////// 


/ 


/,    dU        d(Mu) 


*lw(T*ld 


/////////A 

'< 

77777777777 


An  energy  balance  on  the  gas  and  the  capacitance  combined 
yields , 

dM  =  d(Mu)      ^£ 

ide    de    cde   q~  * 

By  definition  of  the  heat  transfer  resistance  between  the 
capacitance  and  the  surrounding  medium,  R  , 


Cfoo     = 


T    -T  T-T 

Q  00  00 

~~R  =      ~R 
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Now  since  wQ  =  dM/d6  the  expression 

dT    dT  dM      dT 
d6  "  dM  d0  "  W0dM 

can  be  written.   Using  this  expression  and  the  perfect  eras 

relationships  in  the  original  energy  balance, 

T-T 
dM         dM  dT  dT  °° 

-c   T,5tt  +    §?c   T   +   McJs-+c5f+    -5 =    0       • 

p  ld0    d6  v      vd9     cde     R 

c  oo 

Dividing  through  by  c  dM/d9  or  c  wQ  and  noting  that  dT  /d0 

=  dT/de, 

C         T-T 
c  dT         oo 

-kTn  +  T  +  (M+— )§i-  +  5 =  0   . 

1  c   dM  C  R  w~ 

v       v  oo  o 

In  order  to  express  this  equation  in  terms  of  non-dimensional 

temperature,  T*,  it  is  divided  by  T-  to  obtain, 

dT*       dT* 
M*^W  +  C*^tt  -  kT*  +  T*(T*-T*)NTU   =  0   . 
dM*     OdM*      1  °°     0° 

The  variables  are  now  separated  and  the  expression  integrated 

from  the  initial  conditions  (T  =  T~,  M  =  M~)  to  a  state  at  a 

later  time  (T  =  T,  M  =  M) . 

T  M 

T,  fto 

dT* I      dM* 

T*(1+NTU  )-kT*-NTU  T*  =   \  w   M*+C* 
Q  oo     1     °°  °°     \  MQ      0 

T-  '   W 

0  M0 


The  solution  to  the  above  gives , 


(M*+C*) 


log[T*(l+NTUoo)-kT*-NTUT*]    =    -(1+NTUj  lo<?/1+c*\ 
+    log[-kT*-NTU   T*+(1+NTU    )] 

■S    *  00       00  00 
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or 


T*(1+NTU    )     -    kT*    -    NTU   T* 

oo'  2.  OO     00 


(M*+C*) 


(l+c*) 


-(1+NTU    ) 

CO 

[-kT*-NTU   T*+ (1+NTU)] 

CO      oo  J 


Solving   this    expression    for  T*   gives 


fC*+l 


L1+NTU 


T*    = 


(1+NTU    -kT*-NTU   T    )     C*+M*j     +    kT*    +   NTU   T* 

1  oooo'\0/  1  ooco 


1+NTU 

00 

This  expression  may  be  combined  with  the  perfect  gas  law  in 
order  to  obtain  an  expression  for  P*  as  a  function  of  M*. 

.  1+NTU 

/  C  *  +  1  \         oo 


p*  =  m* 


(1+NTU  -kT*-NTU  T  )  C*+M*l   +  kT*+NTU  T* 

oo_L  oo     oo'     I      0  /  0000 

1+NTU 
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Table  1 

Physical  Dimensions  of  Experimental  Apparatus 

1.  Test   Receiver 

Volume:   1.0  3  cubic  feet 
Internal  Area:   5.61  square  feet 
Nominal  Internal  Diameter:   10.19  inches 
Nominal  Outside  Diameter:   10.75  inches 
Nominal  Internal  Height:   21.75  inches 
Mass  of  the  Receiver:   38.5  lbms . 
Material:   6061-T6  Aluminum 
Specific  Heat:   .226  Btu/lbm.  °F 
Thermal  Conductivity:   126  Btu/hr.  ft.  °F 

2.  Flow  Metering  Devices 

3/8  in.  knife-edge  orifice  -  discharge  coeff.  =  .620 
1/8  in.  nozzle  -  discharge  coeff.  =  .858 
1/32  in.  nozzle  -  discharge  coeff.  =  .856 
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Table  2 

Summary 

of  High  Flow 

Rate  Runs 

C*  =  657 

k  =  1.4 

=  1.00 

Run 

Total  Run 

Mass  Flow 

Average 

Average  %  P* 

No. 

Time 
(sec) 

Rate 
(lbm/sec) 

.327 

NTU 

Deviation 

1 

1 

.123 

7.93 

2 

1 

.325 

.092 

8.11 

3 

1 

.325 

.104 

5.59 

4 

1 

.310 

.101 

7.47 

5 

1 

.309 

.089 

6.93 

6 

1 

.309 

.085 

6.11 

7 

1 

.294 

.121 

10.06 

8 

1 

.293 

.084 

5.58 

9 

1 

.293 

.084 

4.92 

10 

1 

.278 

.102 

6.93 

11 

1 

.277 

.105 

7.35 

12 

1 

.277 

.115 

9.63 
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Table  3 

Summary 

of  Low  Flow 

Rate  Runs 

C*  =  657 

Run 

Total  Run 

Mass  Flow 

Average 

Average  %  P* 

No. 

Time 
(sec) 

Rate 
(lbm/sec) 

.0031 

NTU 

Deviation4" 

1 

105 

3.99 

-2.61 

2 

105 

.0031 

2.86 

-3.41 

3 

105 

.0031 

2.90 

-3.18 

4 

105 

.0029 

3.69 

-2.84 

5 

105 

.0029 

3.01 

-3.11 

6 

105 

.0029 

4.78 

-2.71 

7 

105 

.0028 

1.86 

-3.68 

8 

105 

.0028 

3.49 

-2.93 

9 

105 

.0028 

2.18 

-3.39 

10 

105 

.0024 

2.92 

-2.76 

11 

105 

.0025 

3.39 

-2.65 

12 

105 

.0025 

5.54 

-2.22 

13 

105 

.0022 

5.34 

-2.50 

14 

105 

.0022 

7.78 

-1.97 

15 

105 

.0022 

5.61 

-2.09 

16 

105 

.0018 

4.52 

-2.02 

17 

105 

.0018 

7.80 

-1.67 

18 

105 

.0019 

4.43 

-2.01 

19 

105 

.0015 

5.17 

-1.70 

20 

105 

.0015 

8.00 

-1.44 

21 

105 

.0015 

7.91 

-1.27 

22 

105 

.0012 

4.89 

-1.54 

23 

105 

.0012 

6.02 

-1.33 

24 

105 

.0012 

3.83 

-1.70 

25 

105 

.0009 

4.32 

-1.43 

26 

105 

.0009 

6.01 

-1.12 

27 

105 

.0009 

6.21 

-1.07 

28 

105 

.0007 

4.58 

-1.49 

29 

105 

.0007 

3.77 

-1.36 

30 

105 

.0007 

3.47 

-1.45 

31 

90 

.0006 

4.45 

-1.07 

32 

90 

.0006 

13.88 

-0.46 

33 

90 

.0006 

11.19 

-0.52 

34 

75 

.0004 

14.49 

-0.54 

35 

75 

.0004 

14.62 

-0.54 

36 

75 

.0004 

12.39 

-0.51 

A  negative  deviation  indicates  that  the  P*  values  predicted 
by  Reynolds  were  lower  than  those  actually  observed. 
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Table  4 
Summary  of  Isothermal  Sink  Runs 


Run 

Total  Run 

Mass  Flow 

Average 

Average  %  P* 

No. 

Time 
(sec) 

Rate 
(Ibm/sec) 

.0031 

NTU 

Deviation* 

1 

150 

2.69 

-10.39 

2 

150 

.0031 

2.24 

-  8.84 

3 

150 

.0031 

2.17 

-11.08 

4 

150 

.0028 

2.07 

-11.59 

5 

150 

.0028 

2.41 

-10.35 

6 

150 

.0028 

1.97 

-11.47 

7 

150 

.0025 

3.02 

-  9.10 

8 

150 

.0025 

3.30 

-  8.58 

9 

150 

.0026 

2.94 

-  8.94 

10 

9 

.0488 

1.18 

-13.32 

11 

9 

.0478 

1.41 

-13.40 

12 

9 

.0486 

1.38 

-13.22 

13 

9 

.0420 

1.39 

-12.67 

14 

9 

.0422 

1.28 

-13.37 

15 

9 

.0422 

1.15 

-14.34 

16 

9 

.0350 

1.41 

-12.49 

17 

9 

.0348 

1.03 

-13.68 

18 

9 

.0349 

.94 

-13.57 

19 

9 

.0274 

1.73 

-10.78 

20 

9 

.0272 

1.13 

-12.36 

21 

9 

.0272 

1.20 

-12.47 

22 

9 

.0194 

1.85 

-  9.92 

23 

8 

.0192 

1.09 

-11.14 

24 

8 

.0190 

1.05 

-11.14 

25 

7 

.0111 

1.29 

-  7.68 

26 

7 

.0111 

1.20 

-  7.85 

21 

7 

.0114 

1.15 

-  7.98 

A  negative  deviation  indicates  that  the  P*  values  predicted 
by  Reynolds  were  lower  than  those  actually  observed. 


79 


SCHEMATIC  OF   TEST   RECEIVER 


to 

thermocouple 
recorders 


c^n 


-X 


V 

water  level 


large 
/container 


'1 


kto  air 
supply 


/ 


/ 


to  pressure 
transducer 


E 


flow  metering 
device 
mounting 


I  thermocouples 


21  3/4 


10  3/4 


Note:  bleed  valve  not  shown  in  this 

view 


Pig.  1 


80 


Oi 

*— 

a  u 

0)    3    0) 

—  o  -o 

MU    L 

coo 

—  E   u 

in  u   a> 

0)    »- 

.c 

r  *j" 

~- 

- 

C    <U 
C  "O 

-c  O 
o  o 

CM    U 


0) 

> 

> 
C 


O 
10 

J* 

O 

3 
cr 


C 
<0 

a 


u 

4-* 

a 
a 

3 
O 

u 
o 

E 

i_ 

4) 


S 


^ 


0)  -  - 

«/>  — 

a)  a> 

*-»  o 

a) 

L. 

J2 


0) 

>        r 
to 

> 

ai 

0) 


o 

I— 

< 


® 


c 
l_ 


1_ 
0) 

l_Ll 

o 

1 

3 
T3 

o 

c 

CO 

«J 

CM 

u 

4-> 

Z 

• 

0) 

LU 

u. 

3 

t— 

in 

CO 

0) 

> 

a 

CO 

CM 

en 

LU 

c 

(0 


a 
a 

3 


O 
</> 

a> 

V. 

a 

E 
O 

u 


81 


82 


TWiSI 

*\    1 

\M  g§§|            ■ 

« 


\:JJ 


O 


'I-     -  * 


■ 


83 


7.0 


6.5 


6.0 


oS 


5.5 


CO 


cc 


R    5.0 


CO 

to 

<- 


'<■     4.5 


c. 


4.0 


3.5 


/. 

y       Run   1 
+       Run  4 
°       Run  10 

/     + 
4/         n 

+ 

x 

D 

A 

/        C 
/      X 

I- 

/  + 

/      C 

/ 

/  X 
/  + 

□ 

/x 

/+ 

/  n 

3.0 


2.0 


3.0 


4.0 


5.0 


M* 


Fin.    5 
84 


5.0 


4.5 


4.0 


3.5 


in 
■z. 

a: 

c 

_j 

U- 
00 


3.0 


6    2.5 


* 


2.0 


1.5 


1.0 


x/ 

X     Run  3 

+     Run   10 

D     Run   15             > 

a 

+/ 

/ 

/A 

x/ 

V 

! 

i 

2.0  3.0  4.0 

M* 
Fig.    6 


5.0 


85 


7.0 


6.0 


5.0 


CO 


c: 

U. 


to 


4.0 


3.0 


2.0 


1.0 


X 


2.0  3.0  4.0  5.0 

M* 
Fig.     7 


86 


11  .0 


10.0 


* 


U. 
t 

c. 
u 
c 


9.0 


8.0 


7.0 


6.0 


5.0 


4.0 


X 

X        ...  . 

X 

X 

X 

X 

X 

-  

X 

X 

0.09      0.10      0.11       0.12 

NTU  HIGH  MASS  FLOW  RUNS 
Fig.  3 


87 


3.0 


6.0       9.0 
NTU  LOW  MASS  FLOW  RUT'S 


12.0 


Fig.  9 


88 


-15 


-14 


-13 


-12 


U 

c 


c 

L 

c. 


-1 1 


-10 


-9 


-8 


-7 


X 
X 

- 

< 

xx       : 

1 

< 

1 

i 

i 

1 

i 

XX 

X 
X 

X 

1 

I 
1 
1 

X 

X 

X 

rx 

X 

v* 

X 

X 

X 

X 
X 

1.4 


1.8       2.2       2.6 
NTU  ISOTHERMAL  SINK  RUMS 

Fig.  10 


3.0 


89 


10 


NTT 


10 


(D 


10 


-I I I I I  I  t 


10 


8 


10' 

GRPR 
HIGH  MASS  PLOtf  RUNS  1-3 


(1)  laminar  free  convection  line 

(2)  turbulent  free  convection  line 


Pig.  11 


10 


10 

1  =  + 

2  =  x 

3  =  * 


90 


103 


NU 


10' 


(D 


10 


I    I I I— JL 


^L 


1 08  1 09 

GRPR 
HI3H  MASS  FLOW  RUNS  4-6 

(1)  la.ninar  free  convection  line 

(2)  turbulent  free  convection  line 

Pig.  12 


10 


10 

4  =  + 

5  =  x 

6  =  * 


91 


10 


NU 


10' 


(1) 


x 

.IF* 


(2) 


10 


i  i  i  i 


± 


10 


8  109 

GRPR 
HIGH  MASS  PLOrf  RUNS  7-9 

(1)  laminar  free  convection  line 

(2)  turbulent  free  convection  line 

Fig.  13 
92 


10 


10 

7  =  + 

8  =  x 

9  =  * 


10- 


NU 


10' 


(D 


10 


I   I I 111 


10 


8 


10* 

GRPR 
HIGH  MASS  PLOW  RUNS  10-12 


10 


10 


10  =  + 

11  =  x 

12  =  * 


(1 )  laminar  free  convection  line 

(2)  turbulent  free  convection  line 

Pi?.  14 


(2) 


93 


10" 


10 


8 


10* 

GRPR 
LOtf  MASS  PLOW  RUNS  1-3 


1 )  laminar  free  convection  line 

2)  turbulent  free  convection  line 


(2) 


10 


10 

1  =  + 

2  =  x 

3  =  * 


Pig.  15 


94 


10 


10 


10- 

GRPR 
LOtf  MASS  FLOtf  RUNS  4-6 


(1)  laminar  free  convection  line 

(2)  turbulent  free  convection  line 


10 


4  =  + 

5  =  x 

6  =  * 


Fig.  16 


95 


10 


NU 


10' 


(1) 


(2) 


10 


-i i i i  i i  i  i 


JL 


108  109 

GRPR 
LOW  MASS  FLOW  RUNS  10-12 

(1)  laminar  free  convection  line 

(2)  turbulent  free  convection  line 

Pig.  17 


10 


10 
10 

11 

12 


+ 
x 


96 


8 


10w  10v 

GRPR 
LOfl  MASS  FLOtf  HUN 5    15-15 

(1)  laminar  free   convection  line 

(2)  turbulent   free   convection  line 

Pig.    18 


10 


10 


13 

~ 

+ 

14 

= 

X 

15 

~ 

■a 

97 


10- 


10" 

3RPR 
ISOTHERMAL  SINK  RUNS  4-6 


(1)  laminar  free  convection  line 

(2)  turbulent  free  convection  line 


(2) 


4  =  + 

5  =  x 

6  =  * 


Pig.  19 


98 


108  109 

GRPR 

ISOTHERMAL  SINK  RUNS  10-12 

(1)  laminar  free  convection  line 

(2)  turbulent  free  convection  line 

Pig.  20 


(2) 


10 


10 


10 

= 

+ 

11 

= 

X 

12 

~ 

* 

99 


(2) 


GRPR 
ISOTHERMAL  SINK  RUNS  16-1 8 

1 )  laminar  free  convection  line 

2)  turbulent  free  convection  line 

Pig.  21 


16 

~ 

+ 

17 

~ 

X 

18 

~ 

# 

100 


LIST  OF  REFERENCES 

1.  Stanford  Technical  Report  T-l ,  An  Analytical  Investigation 

of  the  Blowdown  and  Charging  Process  in  a  Single  Gas 
Receiver ,  by  W.  C.  Reynolds,  October  1,  1955. 

2.  Stanford  Technical  Report  T-2 ,  An  Experimental  Investiga- 

tion of  Blowdown  Process  in  a  Single  Gas  Receiver, 
Including  the  Effects  of  Heat  Transfer,  by  W.  C. 
Reynolds  ,  J.  W^  Millard,  and  W.  M~i  Kavs ,  December  1, 
1956. 

3.  Grace,  H.  P.  and  Lapple,  C.  E.,  "Discharge  Coefficients 

of  Small  Diameter  Orifice  and  Flow  Nozzles"  Trans- 
actions of  the  ASME,  v.  73,  p.  639-647,  July,  1951. 

4.  Carslaw,  H.  S.  and  Jaeger,  J.  D.,  Conduction  of  Heat  in 

Solids ,  2nd  ed. ,  70-73,  Clarendon  Press,  1959. 

5.  Lyons,  J.  T.,  Pestorius ,  T.  D.,  Sager,  D.  A.,  and  von 

Suskil,  J.  D. ,  Heat  Transfer  in  Closed  Containers  with 
Gas  Injection,  Depart  of  Mech.  Eng.  NPGS ,  Monterey, 
California,  December,  1968. 

6.  Ulrich,  R.  D.,  Wirtz,  D.  P.,  and  Nunn,  R.  H. ,  Thermal 

Analysis  of  Gas  Pressurized  Expulsion  System:  Part  I : 
Heat  Transfer  in  Closed  Containers  with  Gas  Injection, 
Naval  Weapons  Center,  China  Lake,  California,  1968. 


101 


INITIAL  DISTRIBUTION  LIST 


No.  Conies 


1.  Defense  Documentation  Center  20 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

3.  Navy  Ship  Systems  Command  Code  (2052)  1 
Navy  Department 

Washington,  D.  C.   20360 

4.  Asst.  Professor  R.  H.  Nunn,  Code  59Nn  1 
Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Department  of  Mechanical  Engineering  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

6.  Dr.  Richard  D.  Ulrich  1 
Department  of  Mechanical  Engineering 

Brigham  Young  University 
Provo,  Utah   84601 

7.  LTJG  John  T.  Lyons  III  USN  3 
c/o  Mrs.  J.  T.  Lyons 

1991  Grand  Avenue 

St.  Paul,  Minnesota   55105 


102 


TTnnlassi  fieri 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 


(Security  classification  of  title,   body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 

2a.   REPORT    SECURITY    CLASSIFICATION 


i     originating   activity  (Corporate  author) 


Naval  Postgraduate  School 
Monterey,  California   93940 


Unclassified 


2b.  GROUP 


3   REPORT  TITLE 


Heat  Transfer  Considerations  in  a  Pressure  Vessel  Being 
Charged 


4.   descriptive  NO  T  E  s  (Type  of  report  and,  inclusive  dates) 

Master's  Thesis;  (June  1969) 


3     AUTHOR(S)  (First  name,  middle  initial,  last  name) 


John  Thomas  Lyons  III 


6.    REPORT    DATE 


June    1969 


7a.    TOTAL    NO.    OF   PAGES 


102 


7b.    NO.    OF    REFS 


8«.    CONTRACT    OR    GRANT   NO. 


b.    PROJECT   NO. 


9a.    ORIGINATOR'S    REPORT   NUMBER(S) 


9b.  OTHER   REPORT  noisi  (Any  other  numbers  that  may  be  assigned 
this  report) 


10.    DISTRIBUTION    STATEMENT 


Distribution  of  this  document  is  unlimited, 


II.    SUPPLEMENTARY    NOTES 


12.    SPONSORING   MILITARY    ACTIVITY 


Naval  Postgraduate  School 
Monterey,  California   93940 


13.  ABSTRACT 


Experimental  data  for  the  charging  of  an  air  receiver  is 
presented  and  interpreted  in  detail.   The  data  indicates  a 
substantial  departure  from  the  adiabatic  behavior.   The 
experimental  results  are  used  to  evaluate  existing  closed 
form  expressions  for  the  thermodynamic  state  of  a  gas  in  a 
receiver.   A  method  for  experimentally  determining  the  con- 
vective  heat  transfer  coefficient  is  developed,  evaluated  and 
used  in  conjunction  with  these  expressions. 

The  experimental  work  was  performed  from  March  19  69 
through  May  1969  at  the  Naval  Postgraduate  School,  Monterey, 
California. 


DD 


"""*  1473 

i  nov  es  I  "T  /  sJ 
S/N    0101-807-681 1 


(PAGE    1) 


103 


Unclassified 


Security  Classification 


A-3140* 


Ilnclass.     'ip.d 


Security  Classification 


KEY     WORDS 


Convective   Heat  Transfer 
Charging  of  Receiver 
Thermodynamic   State   of   a  Gas 


DD 


FORM 

I    NOV   69 


1473  <BACK> 


Unclassified 


S/N     010I-907-60?! 


104 


Security  Classification 


thesL9947   ) 

Heat  transfer  considerations  in  a  pressu 


3  2768  002  12451  3 

DUDLEY  KNOX  LIBRARY 


